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ABSTRACT 

We present a survey of the extinction properties of ten lensing galaxies, in the redshift range z = 
0.04-1.01, using multiply lensed quasars imaged with the ESO VLT in the optical and near infrared. 
The multiple images act as 'standard light sources' shining through different parts of the lensing galaxy, 
allowing for extinction studies by comparison of pairs of images. We explore the effects of systematics 
in the extinction curve analysis, including extinction along both lines of sight and microlensing, using 
theoretical analysis and simulations. In the sample, we see variation in both the amount and type of 
extinction. Of the ten systems, seven are consistent with extinction along at least one line of sight. The 
mean differential extinction for the most extinguished image pair for each lens is A{V) = 0.56 ± 0.04, 
using Galactic extinction law parametrization. The corresponding mean Ry = 2.8 ± 0.4 is consistent 
with that of the Milky Way at Ry = 3.1, where Ry = A(y)/E{B - V). We do not see any strong 
evidence for evolution of extinction properties with redshift. Of the ten systems, B1152+199 shows the 
strongest extinction signal of A{V) = 2.43 ±0.09 and is consistent with a Galactic extinction law with 
Rv — 2.1 ± 0.1. Given the similar redshift distribution of SN la hosts and lensing galaxies, a large 
space based study of multiply imaged quasars would be a useful complement to future dark energy 
SN la surveys, providing independent constraints on the statistical extinction properties of galaxies 
up to z 1. 

Subject headings: dust, extinction — galaxies: ISM — gravitational lensing 



1. INTRODUCTION 

The study of extinction curves of galaxies at high 
redshift has ge nerated a l ot of i nterest in recent 
years (see e.g.. IRiess et ahl l199fia[ iFaIco et l1999t 
iGoudfrooii' '2000* 'Muruhv & Liske' '20041; 'Kann et all 
11)06: Goicoechca ct al. 2005; York et al. 2006). Light 
reaching us from distant sources is extinguished by dust 
along its path making it important to correct mea- 
surements for the amount and properties of the ex- 
tinction. Extragalactic dust extinction can for exam- 
ple affect measurements of Type la supernovae (SNe 
la) u sed to determine various cosmological p arameters 
(e.g.. IRiess et al.lll998l IPerlmutter et alJ 11999^1 and the 
star-formation rates for high redshift starburst galaxies 
which are used as probes of galaxy evolution (see e.g., 
iMadau et"a l. 1998). Yet, even if dust properties and thus 
extinction may vary with redshift and environment, an 
average Galactic extinction law is often applied when cal- 
ibrating extragalactic data due to the lack of knowledge 
of the extinction properties of higher redshift galaxies. 

Traditionally, extinction curves are measured by com- 
paring the spectra of two stars of the same spectral type, 
which have been re ddened by different amounts (see e.g., 
iMassa et al.l 119831) . As it becomes significantly harder 
to measure spectra of individual stars with distance 
this method is limited in its application to the Milky 
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Way and the nearest galaxies. The extinction curves of 
the Milky Way along different lines of sight have been 
mapped extensively using this method and have been 
shown to follow an empirical parametric function which 
depends only on one parameter, Ry = A{V)/E{B — V), 
where A{X) is the total exti nction at wavelength A and 
E{B -V) ^ A{B) - A(y) iCardelh et & V^19m . The 
mean value of Ry in the Milky Way is 3.1 fCardell i et alJ 
119891) but for different lines of sight the value rang es from 
as low as Ry « 1.8 toward the Gal actic bulge (lUdalsk) 
2 90^ and a,s high as K ,. « 5.6-5.8 llCardelli et 311119891: 
Fitzpatrick et al.lll999D . A lower Ry corresponds to a 
steeper rise of the extinction curve into the UV, whereas 
it has little effect on the extinction in the infrared. 

Extinction curves have also been obtained for the Small 
and Large Magellanic Clouds (hereafter, SMC and LMC, 
respectively) and M31 using this method. The mean ex- 
tinction curve of the LMC differs from the Galactic ex- 
tinction law in that the bump at 2175 A is smaller by 
a factor of two (as measured by the residual depth of 
the bump when the continuum has been extracted) and 
the curve has a steep er rise into the UV for wavelengths 
shorter than 2200 A l|Nandv et alJll98lD . The extinction 
curve of the SMC is well fitted by an A{X) cx A^^ curve 
which deviates significantly from the Galactic extinc- 
tion law and the LMC extinction for A^^ > 4 /im~^ and 
in particular s hows no bump at 2 17 5 A and a steep er rise 
into the UV l|Prevot et al.lll984|) . iBianchi et alT l|1996D 
found that the extinction of M31 follows that of the av- 
erage Galactic extinction law. The various extinction 
properties shown by these galaxies, especially in the UV 
and shorter wavelengths, further strengthens the need to 
find a method to study the extinction curves of more 
distant galaxies. 

A few methods have been proposed for measuring ex- 
tinction curves for more distant galaxies. One method 
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is basically an extension of the traditional method of 
comparing stars of the s ame spectral type to coniparin K 
the SNe la ijRiess et al.lll996at IPerimutter et all 119971) . 
The extinction is estimated from comparison with unred- 
dened, photometrically similar SNe la. A subset of SNe 
la, with accurately determined extinctions and relative 
distances, is then used to further determine the relation- 
ship between light and color curve shape and luminos- 
ity in the full sample. SN la extinction studies usu- 
ally give lowe r Ry values than t he mean Galactic value 
of Rv = 3.1 llRiess c t al. 1996bt iKrisciunas e'ram2000l 
IWang 200^). 

Quasars with damped Lya syster ns fPLAs] in the fore- 
ground have also been studied bv iPei et alJ l|199H) and 
were found to be on average redder than those without. 
By comparing the optical depths derived from the spec- 
tral indices and the ones derived from excess extinction 
at the location of the Galactic extinction law bump they 
found that their sample of five quasars with DLAs is not 
consistent with the Galactic extinction law, marginally 
compatible with the LMC extinction and fu lly compati- 
ble with SMC extinction. iMurphv fc Liskd 1(20041) stud- 
ied a larger sample of the Sloan Digital Sky Survey 
quasars with damped Lya systems in the foreground 
and found no sign of extinction. They suggest that the 
differe nce between their results and that of Pci et alJ 
may be due to the s r nall number statisti cs in 
the study bv lPei et alJ l|199H) . lEllison et all l|2005|) also 
found that intervening galaxies cause a minimal redden- 
ing of background quasars in agreement with the re- 
sults of Murohv & Liskc ( 2004) while York et al. ( 2006) 
found E{B — V) of up to 0.085 for quasars from the 
Sloan Digital Sky Survey with Mg II absorption. In their 
studv lYork et all (||006) fou nd no evidence of the 2175 A 
bump (at variance with iMalhotra I (|1997|) ') and found 
that the extin ction curves are similar to SMC extinction. 
lOstman et al.l |2006) studied the feasibility of measuring 
extinction curves by using quasars shining through galax- 
ies. For the two such systems which survived their cuts, 
they argued that the extinction curves in the foreground 
spiral galaxies were consistent with Galactic extinction. 
They further suggested a possible evolution in the dust 
properties with redshift, with higher z giving lower Ry 
by studying values obtained from the literature in addi- 
tion to their own. 

Extinction curves of high redshift galaxies have also 
been studied by looking at the spectral energy distri- 
bution of gamma-ray bursts (GRBs). For example, 
ITakobsson et al.l (|2004|) fitted a Galactic extinction law, 
an SMC and an LMC extinction law to the afterglow of 
GRB 030429. The afterglow, at z = 2.66, was best fit by 
an SMC lik e extin ction curve with A(y) = 0.34 ± 0.04. 
iKann et alJ (|2006|) studied the extinction of a sample of 
19 GRB afterglows and fitted them to various dust ex- 
tinction models. They found that the SMC extinction 
law was preferred by a great majority of their Golden 
Sample (seven out of eight) while one afterglow was 
best fit by a Galactic extinction law (the other eleven 
were equally well fit by SMC, LMC and Galactic ex- 
tinction). The mean extinction in the y-band was 
A (V) = 0.21 d z 0.04 . 

Goudfrooii (2000) reviewed the dust properties of giant 
elliptical galaxies and found that they are typically char- 
acterized by small Ry if they are in the field or in loose 



groups, but that if they are in dense groups or clusters 
their Ry values are close to the mean Galactic Ry = 3.1. 
Earlytype elliptical galaxies typically have low A{V) (see 
e.g.,Epudfrooii 199i, who found A{V) < 0.35 for dust 
l anes in elliptica l s). 

iNadeau et al.l (|199iri pointed out that gravitationally 
lensed quasars could be used to mea sure the ext i nction 
curves of higher redshift galaxies. iFalco et all (|1999D 
explored a large sample of 23 lensing galaxies using 
this method and found that only seven were consistent 
with no extinction. This meth od has also b e en ap - 
plied to single sy s tems by e.g. I.Taunsen et al.l ll l997t): 
iToft et abl mm-. IMotta et alJ J2002I): IWucknitz et aJJ 
(120031): iMufioz et alJ J2004D : iWisotzki et alJ (|2004|) : 
lGoicoeche^^r^.1 (|2005D and shows varying extinction 
properties between different lensing systems. 

Here we present a systematic study of the extinction 
curves of gravitational lenses based on a survey of 10 
lens systems. We have made a dedicated effort to mini- 
mize the number of unknowns and effects that can mimic 
extinction. We have broad wavelength coverage in nine 
different optical and NIR broad bands. An effort was 
made to minimize the time between the observations for 
each system in the different bands to minimize the ef- 
fect of intrinsic quasar variability and microlensing. All 
our systems have spectroscopically determined redshifts 
for both the quasar and the lensing galaxy. Finally, our 
systems span the range oi z = 0.04-1.01 giving us the 
possibility to explore possible evolution with redshift. 

The outline of this paper is as follows: In §[3 we de- 
scribe the details of the employed method and discuss dif- 
ferent sources of systematics and of random errors which 
may affect our results. We also present the results of 
simulations which explore the effects of these errors on 
data sets similar to those we obtain in the survey. In § |21 
we present the data and the data reduction of the ESQ 
VLT survey for the 10 lensing systems. We present the 
results of our analysis of each individual system in § 14.11 
and the analysis of the full sample in § 14.21 Finally we 
summarize our results in § 

2. METHOD AND SIMULATIONS 

In this section we introduce the method and explore 
the different sources of systematic and random errors 
through simulations. In particular, we study the effects 
of achromatic microlensing and of extinction along both 
lines of sight, and use simulations to explore the con- 
ditions, under which it will be possible to recover and 
distinguish between different extinction laws. 

2.1. Lensing 

Gravitational lensing is the deflection of light rays 
due to the gravitational field of the matter distribution 
through which the light passes. For a geometrically thin 
lens, i.e., where the depth of the lens is small compared 
to the distance between the lens and observer, Di, and 
the lens and the light source. Dig, the deflection angle is 
given by 

a(r) = lj /dVl](r)^^ (1) 
7 |r — r p 

where G is the gravitational constant, c is the speed of 
light, I](r) is the surface mass density of the lensing mass 
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Fig. 1. — The lensing setup. D; is the distance to the lens, Ds 
the distance to the source, Di^ the distance from the lens to the 
source, a is the deflection angle (see eq. s is the position vector 
of the source in the source plane, and r is the position vector in the 
lens plane. The figure shows a doubly imaged source, where one 
line of sight goes through the lensing galaxy and the other passes 
outside it. 



and r is the impact vector of the hght ray (see Figure ^ 
for a sketch of the lensing setup). For a position vector 
s in the source plane one will see images at locations r 
in the lens plane which satisfy the lens equation 



D, 



r - Dis&ir) 



(2) 



where Dis is the angular diameter distance between the 
lens and the source and where we have assumed that 
the size of the lens is small compared to Di , Dg and Dis . 
When eq. (0) has more than one solution we see multiple 
images of the source in the lens plane which are in general 
located at different distances from the center of the mass 
distribution. In the case of distant quasars being lensed 
by foreground galaxies, the condition that the lens size 
be small compared to Di, Ds and Dis is fulfilled. 

The multiple images of lensed quasars act as 'standard 
candles' shining through different parts of the lensing 
galaxy and can therefore be us ed to studv its extin ction 
curve (as first pointed out bv lNadeau et al.lll99H) . As 
lensing is achromatic in nature, the flux ratio of any two 
lensed images should be independent of wavelength. If, 
however, one of the images shines through a dusty part 
of the galaxy and the other image does not, the first 
image will appear red compared to the other. By map- 
ping the flux ratio as a function of wavelength one can in 
principle directly trace the differential extinction curve 
between the two images, without making any assump- 
tions about the intrinsic spectrum of the quasar (as it 
cancels in the calculation of the flux ratios). Depending 
on the number of images, it will be possible to obtain dif- 
ferential extinction curves for several paths through the 
lensing galaxy. 

As the light rays travel along different paths for mul- 
tiple images, their travel times do in general differ, in- 



troducing a time delay between different images. If the 
quasar is variable, the variability will show up at different 
times in the different images which can lead to inaccurate 
estimates of the dust extinction. Ideally one would like 
to measure each image with a time separation according 
to the time delay to correct for this effect. Time delays 
are, however, difficult to measure and measurements ex- 
ist for only a few lenses. Alternatively one can observe 
simultaneously in all the observing bands. This would 
mean that any achromatic variability would cancel out 
when comparing the images (but would lead to biased 
estimates of the intrinsic brightness ratio of the images). 
Simultaneous observations also have the additional ben- 
efit that the effects of achromatic microlensing will, to 
first order, only affect the intrinsic ratio estimate and 
not the shape of the extinction curve, which would in 
general not be the case if the images were observed ac- 
cording to a time schedule given by the time delay. The 
effects of microlensing are the greatest potential source 
of systematic error in our extinction curve analysis and 
are addressed in greater detail in § 12.41 

2.2. Extinction 

As lensing is an achromatic process we would expect 
the magnitude difference in all bands to be constant for 
each pair of lensed images in the absence of extinction. 
Extinction reduces the brightness of the measured images 
by a different amount for each band (and image) depend- 
ing on the amount and properties of the dust along the 
line of sight to the images. It is this difference which gives 
rise to the extinction curve as a function of wavelength. 
As both the images might be affected by extinction what 
one is really measuring is the differential extinction be- 
tween the pair of images. The extinction affects each 
measured data point as: 



m{X) =m{X)+A{X), 



(3) 



where to(A) is the measured magnitude of the image, 
m(A) is the intrinsic magnitude of the image and A(X) is 
the extinction at wavelength A. When comparing images 
A and B one therefore gets: 



m^(A) 



ATO(A) = m^(A) 

= [m^ - m^) + {X) 
= A7h + A'^'ff{X}, 



(4) 



A^{X) 



where Am = — rh^ is the intrinsic magnitude 
difference which does not depend on wavelength and 
^d«//(A) = A^{X) - A^(A) is the effective differential 
extinction law as a function of the wavelength. 

We consider three different extinction laws and assume 
that the extinction of one of the images dominates the 
other (see further discussion on extinction in both im- 
ages in § I2.3|l . The first extinction law we consider is 
the empirical G alactic extinction law as parametrized by 
ICa,rde]li et ^ l(T98^ : 

(5) 



A{X) = E{B - V) [Rva{x) + b{x)] 



:A{V) 



Hv 



where A{X) is the total extinction at wavelength A, 
E{B ~V) ^ A{B) - A{V) is the color excess, Ry = 
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A{V)/E{B — V) is the ratio of total to selective extinc- 
tion, a{x) and b{x) are polynomials and x — X~^. We 
also consider an extinction law which is linear in inverse 
wavelength which is characteristic for the extinction in 
the SMC: 

Finally we extend the linear law to a power law: 

A 



^(A) = A{V) 



5500A 



(7) 



To fit data points to the extinction laws one first shifts 
the wavelength of the measured bands to the rest frame 
of the lensing galaxy, i.e. Xj — A^/(l + zi) where is 
the observed wavelength in band j and z; is the redshift 
of the lensing galaxy. For each image pair one can then 
calculate the magnitude difference of the images in each 
measured band: 



Ha, 



-2.5 log 




(8) 



where / fj^ is the flux ratio between images labeled B 
and A at Xj . One can then perform fits for eq. lO where 
one replaces A'^^ff{V) with A{V): 

m^{X)^m^{X) = A{X) + Am, (9) 

where ^(A) is one of the extinction laws described in 
eqs. ©, © and O. If radio measurements exist for 
the flux density of the images, one can use them to con- 
strain the intrinsic magnitude difference, Am, as radio 
measurements are not affected by extinction. 

2.3. Extinction along both lines of sight 

As the method we use measures differential extinction 
curves, we wish to investigate the systematics of extinc- 
tion along both lines of sight in our results. We therefore 
study the effects of a Galactic extinction law along both 
lines of sight but with different values of Rv- When 
both images suffer extinction we expect to get an effec- 
tive extinction law which may have different properties 
to those of either line of sight. In the general case, when 
E^{B -V) ^ E^{B -V), the difference in extinction 
suffered by image A vs. image B will be given by an 
effective Galactic extinction law: 

A'^'ff{X)=A^{X) - A^{X) (10) 
= E^{R^a{X-') + biX-'))- (11) 
^^(i?^a(A-i) + fe(A-i)) 



E^ 



fS^ - E^) 



Hj JXy — III JXy 



EB -E^ 



a{X-') + b{X-m) 



(E^ - E^ 



) {Rv"aiX- 



b{X- 



where we have written E(B ~ V) ^ E for simplicity, 
explicitly used the assumption E^ ^ E^ in the step 
from equation ^ and ^ and Ry-'^^ = [E^ R^ - 
E^Ry)/{E^ - E^) is the effective Ry we measure. 
For completeness we note that in the special case of 
E^ — E^ the resulting effective extinction curve is not 
given by the Galactic extinction law parametrization. If 




Fig. 2. — Contour plot of r; as defined in eq. 1141 which measures 



the relative error of the effective R' 



■diff 



to 



as a function of 



E^/E^ and Ry/Ry for extinction along both lines of sight. 



we take E^ > E^, we find that the ratio of Ry^^ to the 
W 



Ry (which we wish to measure) is given by: 



R 



dtff 



E^ -E^R^/R 



Ri 



I 



EB -E^ 
E^ 



(13) 



1 



Ry 



EB-E^ 
= 1 + 11 

so the error introduced in our estimate due to the non- 
zero extinction of image A is 



E^jE^ 
l-E^/E^ 



1 - 



Ry 

Ry 



(14) 



We note that if Ry > i?y, the inferred value of Ry for 
image B will be lowered, and vice versa, and that, in 
theory, any value of Ry can be obtained. A contour 
plot of ri can be seen in Figure |21 If we take the most 
extreme values of the Milky Way (2 < Ry < 6) we see 
that the bracket in equation 114|l can realistically range 
from (1 - 6/2) = -2 to (1 - 2/6) = 2/3. In the worst 
case scenario, when Ry ^ i?^, we therefore need the 
ratio of E'^/{E^ - E^) to be half that of the desired 
accuracy, i.e., for a desired accuracy of 10% in Ry we 
need E^/{E^ — E^) < 0.05. For more realistic cases of 
Ry = 3 for image B and Ry = 4 for image A we only 
need E^/{E^ — E"^) < 0.3. In general, for an accuracy 
of \ri\ < 770, we need: 



El 

EB 



< 



Vo/\l 



Ry/Ryl 



l + ?7o/|l 



Ry/Ryl 



(15) 



Finally we note for completeness that a linear extinc- 
tion in both images trivially produces a linear differen- 
tial extinction whereas a power law extinction along both 
lines of sight does not in general produce a power law for 
t he differential extinctio n. 

McGough et alJ l)2005tl did a similar study of the effect 
of non-zero extinction along both lines of sight, where 
they also found that, in theory, any value of Ry can be 
obtained. They suggest that cases where only one of 
the images is lightly reddened or the dust properties are 
the same for both sight lines are likely rare and hard to 
confirm. We point out that it is not crucial that one of 
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the images have no or Uttle extinction in absolute terms, 
but only relative to the image we are comparing it to, 
and that in practice, one of the images often shows less 
extinction than the others. When dealing with multiply 
imaged quasars, in particular for doubly imaged systems, 
one of the images is often situated at a greater distance 
from the lens galaxy than the others and therefore may 
be less affected by extinction. In some cases this lack of 
extinction in one of the components can be confirmed by 
studying the images in the X-rays (K. Pedersen et al., 
2006, in preparation). 

2.4. Microlensing 

Microlensing, lensing by stars or other compact objects 
in the lens galaxy, can also affect our data and in par- 
ticular it can affect the continuum part of the emission. 
This is becau se, according to standard quasar models 
llKrolildll999D . the regions giving rise to the continuum 
and emission lines are of different size and therefore af- 
fected differently by microlensing which is more effective 
on small scales. As the emission lines arise from regions 
several orders of magnitude larger than the region emit- 
ting the continuum, the microlensing acts strongest on 
the continuum emission but should be nearly absent for 
emission lines. 

To allow for possible corrections due to this effect we 
calculate for each quasar the ratio of the spectral line 
emission to the total emission in each band (this varies 
for the systems as the quasars are at different redshifts) . 
We use a composite quasar spectrum for our calculations 
as derived byiVanden Berk et al. (2001.) using 2200 spec- 
tra from the Sloan Digital Sky Survey. In accordance 
with their results we model the continuum as a broken 
power law {fx oc A"^) with ax = -1.56 for A < 4850 A 
and ax = 0.45 for A > 4850 A. We take into account 
spectral lines with equivalent width W > 1 A and model 
them as Gaussians which we add to the continuum to get 
the final template spectrum. Using this standard quasar 
template we calculate the ratio of the flux coming from 
the continuum compared to the total emission for each 
measurement band shifted to the corresponding band at 
the redshift of the lensing galaxy using the transmission 
curves of the corresponding filter (see Figure I^J. 

We add to our fit an effect from an achromatic mi- 
crolensing signal which affects the different bands propor- 
tionally to the ratio of the continuum to the total emis- 
sion (consisting of the continuum with the added spectral 
lines). Achromatic microlensing affects the fluxes of the 
images as fj — > /j(l -I- stj) where rj is the ratio of the 
continuum emission to the total emission in band j and s 
is a microlensing parameter giving the strength of the mi- 
crolensing signal (which is a constant in the achromatic 
case). The measured magnitude is therefore: 

m(A,) = -2.51ogio im + sr,)) + A{X,) (16) 
= m(Aj) - 2.51ogio(l + sr,) + A{X,). 
The additional term in eq. (|16|l modifles eq. ^ to: 



(m^ - m^)(Aj) = A(Aj) + Am - 2.51ogi 



1 



s^r. 



.A(A,) 



2.5 



In(lO) 




Am - 



1(A) 



Fig. 3. — The composite quasar spectra modeled as a broken 
power law with spectral lines with equivalent width W > 1 A 
added on as Gaussians. The superimposed transmission curves 
correspond to the FORSl and ISAAC U, B, V, R, I, z (Gunn) , z 
(special), Js, H and Ks bands. The transmission curves have been 
shifted to show which part of the spectrum they correspond to for 
a quasar of 2 = 2.2. 



where s = 2.5(s"^ — s^)/ln(10) and the approximation 
is made to reduce the number of parameters in the flt. 
For the approximation in eq. H17() to be valid we need 
rj\s^ — < 1. The value of Vj can in theory lie be- 
tween and 1 but in our case lies between 0.8 and 1, so 
therefore we need — s"^| < 1 in practice. In physical 
terms this is roughly equivalent to the condition that the 
change in magnitude difference between the two images 
due to microlensing should be less than 1 mag. Most mi- 
crolensing studies show magnitude changes of less than 
0.5 mag so therefore we do not expect this approximation 
to affect our results. 

Microlensing can also introduce a smooth chromatic 
signal similar to an extinction signal if the source has 
different colors in different locations. Our methods how- 
ever do not take this into account as our broad band 
photometry does not contain enough information to dis- 
entangle such an effect from extinction. To separate the 
two it is necessary to study the spectra of the object it- 
self (and not a mean q uasar spectrum) along the lines of 
iWiicknitz et"aD li200l . 

2.5. Monte Carlo Analysis 

We estimate the errors of the parameters of the fit- 
ted extinction curves using Monte Carlo simulations. We 
generate a thousand realizations of the data set by allow- 
ing each data point to vary with a Gaussian distribution 
centered on the measured point and with a variance cor- 
responding to the error estimate of the data point. We 
then apply the extinction law fits separately to each data 
set. If the fits are good, i.e., there is no degeneracy or 
systematic errors, the resulting fitted parameters should 
also follow a Gaussian distribution. We quote the median 
of the parameter distribution as our 'best fit parameter' 
and the quoted error is the standard deviation of the pa- 
rameter distribution. For the Galactic extinction law we 
constrain < Ry < 7 unless otherwise noted. This 
method is applied to both the simulations in § 12.61 and 
the real data sets in § ^ 

2.6. Simulated Data 

To test the ability of our analysis to recover and dis- 
tinguish between different extinction curves, and the ef- 
fects of microlensing and noise, we apply our method 
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to simulated data. Our simulated data consist of data 
in nine bands, UBVRIzJsHK, with an applied signal 
corresponding to extinction, noise and microlensing. We 
do runs for different kinds of the three extinction laws, 
with or without noise, and with or without microlensing 
effects. 

2.6.1. Pure extinction in one image 

The first set of runs we do are pure extinction in one 
image with a 0.05 mag la uncertainty in each data point. 
The purpose of those runs is twofold, first to test our 
routines and secondly to see whether there is a significant 
difference between the goodness of fit for the different 
extinction laws. We find that in all cases our routines 
converge to the given initial parameters. In addition we 
find that the ability to recognize one extinction law from 
the other depends on the strength of the extinction and 
the redshift of the simulated data. This is not surprising 
as the three fitted extinction laws behave very similarly 
for A > 2500 A and for z = all our bands lie above that 
limit. The importance of the strength of the extinction 
is also easy to understand, as the ability to detect any 
difference between extinction laws will be overwhelmed 
by the photometric errors for very weak extinction. 

At z = for a Galactic extinction law input the re- 
duced chi-squared, xt^ for ^ Galactic extinction law is 
somewhat lower than for the other extinction laws (~ 1.0 
for the Galactic extinction law, w 1.0-1.6 for the power 
law, « 1.0-6 for the linear law). The difference de- 
pends both on the strength of the extinction, A(V) (with 
A{V) < 0.1 resulting in equal goodness of fits), and on 
the value of Ry with more extreme values (Rv — 1, 
Rv = 6) giving a greater difference than Ry = 3.0 (see 
Table n] for representative values and Figures 0] and for 
representative plots). When the input extinction law is a 
linear law, A~^, the power law yields the same goodness 
of fit {xl ~ 1-0) but the Galactic extinction law gives 
a slightly higher value (« 1.0-1.9). Finally, when the 
input extinction law is a power law with power index 
a = 2, there is a significant difference in the goodness 
of fit with the power law fit giving xt ~ l-O but the 
Galactic extinction law giving xt ~ 1-0-1.7 and giv- 
ing xi ~ 1.0-6.5. Here the difference in the goodness of 
fit for the Galactic extinction law and A^^ depends on 
the value of A{V), with lower A{V) (i.e., less extinction) 
giving lower xt for the two extinction laws. 

We also run simulations for two non-zero values of the 
redshift, at z — 0.30 and z = 0.80. At these higher red- 
shifts the difference in the extinction laws becomes more 
prominent which is reflected in the goodness of the fits 
(see Tableland Figs. 0] and [Sjl making them in general 
easier to distinguish from one another. This is because 
the extinction laws are all very similar for A~^ < 3 ^m~^, 
and for redshift of z = our lowest wavelength band, the 
[/-band, lies below this limit. At the higher redshifts, the 
lowest wavelengths move into the UV, which is more sen- 
sitive for differences in the extinction laws. For very weak 
extinction {A{V) < 0.1) it is still the case that the dif- 
ferent types of extinction laws become hard to separate, 
as the error bars can dominate the extinction signal. 

From our simulations we can deduce that given a 
strong enough extinction the different kinds of extinc- 
tion laws should be recognizable from each other. As 
we can see from Tabled the needed strength is depen- 
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Fig. 4. — A sample plot for simulations showing three data sets 
at redshift of 2 = 0.0, 0.3 and 0.8. The extinction law applied is 
the Galactic extinction law with Ry = 3 and A{V) = 0.1, 0.5, 1.0 
and the three fits shown are the Galactic extinction law (dashed 
line, eq. power law (dotted line, eq. J7J) and linear law (solid 

line, jnj). 
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Fig. 5. — A sample plot for simulations showing three data sets 
at redshift of z = 0.0, 0.3 and 0.8. The extinction law applied is 
the power law with a = 2 and AiV) = 0.1,0.5, 1.0 and the three 
fits shown are the Galactic extinction law (dashed line, eq. JFJ), 
power law (dotted line, eq. Q) and linear law (solid line, ISl'l. 
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TABLE 1 
Goodness of fits from simulations 



Input parameters Output xi from the different fits 



z 


Type 


AiV) 


Parameter 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0.0 


Galactic 


n 

u. 


c; 

■ O 


ity 






0.94 


0.93 


1.6 


1.6 


5.9 


3.5 


0.3 


Galactic 


n 


c; 

■ O 


ity 




1 n 

J- .u 


0.97 


0.96 


1.9 


1.6 


8.5 


4.2 


0.8 


Galactic 


n 


c; 

■ O 






1 n 

J- .u 


0.95 


0.94 


3.4 


2.9 


15 


6.8 


0.0 


Galactic 


l 


.0 






3.0 


0.95 


0.94 


1.4 


1.1 


2.5 


1.1 


0.3 


Galactic 




n 


ity 




o.u 


0.96 


0.96 


1.9 


1.2 


2.3 


1.5 


0.8 


Galactic 


1 

i, 


n 
.U 


Rt r 




o.u 


0.96 


0.94 


1.8 


1.7 


2.9 


1.6 


0.0 


Galactic 


u. 


r. 
.0 


Rt r 

JXy 




o.u 


0.95 


0.94 


1.1 


1.0 


1.5 


1.0 


0.3 


Galactic 


n 


c; 

■ O 


R^r 

ixy 




o.u 


0.95 


0.95 


1.2 


1.0 


1.4 


1.1 


0.8 


Galactic 


n 


c; 

■ O 


Rt. 

ity 




o.u 


0.96 


0.94 


1.2 


1.2 


1.6 


1.2 


0.0 


Galactic 


n 




Rtt 

riy 




o.u 


0.97 


0.97 


0.95 


0.94 


0.98 


0.92 


0.3 


Galactic 


n 




Rt. 

ity 




o.u 


0.95 


0.94 


0.97 


0.98 


0.97 


0.96 


0.8 


Galactic 


n 




Rt. 

ity 




o.u 


0.95 


0.95 


0.97 


0.98 


0.99 


0.96 


0.0 


Galactic 




■ O 


Rtt 

Ity 




u.u 


0.96 


0.96 


1.3 


1.0 


1.2 


1.3 


0.3 


Galactic 


n 

u. 


c; 

■ O 


Rt. 

ity 




u.u 


0.96 


0.96 


1.7 


1.3 


1.9 


1.8 


0.8 


Galactic 


n 

u. 


c; 

■ O 


Rtt 

Ity 




u.u 


0.95 


0.94 


1.6 


1.3 


2.7 


1.9 


0.0 


Power law 


l 


.0 


a — 


2 




1.6 


1.7 


0.96 


0.95 


6.5 


3.5 


0.3 


Power law 




.0 


a — 


2 




2.3 


2.1 


0.95 


0.95 


11 


5.7 


0.8 


Power law 




n 


a — 


2 




4.7 


4.5 


0.95 


0.93 


21 


11 


0.0 


Power law 


u. 


c: 
.0 


Ot = 


9 




1.1 


1.2 


0.97 


0.96 


3.4 


1.9 


0.3 


Power law 


n 


■ O 


a — 


2 




1.4 


1.3 


0.96 


0.95 


5.5 


3.0 


0.8 


Power law 


n 

u. 


c; 

■ O 


a — 


2 




2.4 


2.3 


0.94 


0.94 


11 


5.5 


0.0 


Power law 


n 


Q 

.o 


a — 


2 




1.0 


1.0 


0.97 


0.97 


2.2 


1.4 


0.3 


Power law 




Q 
.O 


a — 


2 




1.1 


1.1 


0.96 


0.95 


3.4 


1.9 


0.8 


Power law 


n 


Q 
■ O 




2 




1.6 


1.6 


0.95 


0.94 


6.4 


3.4 


0.0 


Power law 


0, 


,1 


a = 


2 




0.95 


0.95 


0.94 


0.94 


1.2 


1.0 


0.3 


Power law 


0, 


,1 


a = 


2 




0.97 


0.96 


0.95 


0.96 


1.4 


1.1 


0.8 


Power law 


0, 


,1 


a = 


2 




1.0 


1.0 


0.94 


1.0 


2.3 


1.4 


0.0 


Linear 


1, 


,0 


a = 


1 




1.9 


1.1 


0.96 


0.94 


0.96 


0.96 


0.3 


Linear 


1, 


,0 


a = 


1 




2.2 


1.2 


0.95 


0.95 


0.96 


0.95 


0.8 


Linear 


1, 


,0 


a = 


1 




1.9 


1.5 


0.95 


0.95 


0.96 


0.95 


0.0 


Linear 


0, 


,5 


a = 


1 




1.3 


1.0 


0.96 


0.96 


0.96 


0.96 


0.3 


Linear 


0, 


,5 


= 


1 




1.4 


1.0 


0.96 


0.94 


0.96 


0.95 


0.8 


Linear 


0, 


,5 


ct = 


1 




1.3 


1.1 


0.94 


0.93 


0.96 


0.94 


0.0 


Linear 


0, 


,3 


a = 


1 




1.1 


0.96 


0.96 


0.95 


0.96 


0.96 


0.3 


Linear 


0, 


,3 


CK = 


1 




1.1 


0.97 


0.96 


0.95 


0.95 


0.95 


0.8 


Linear 


0, 


,3 


ct = 


1 




1.1 


1.0 


0.94 


0.92 


0.95 


0.94 


0.0 


Linear 


0, 


,1 


CK = 


1 




0.98 


0.95 


0.95 


0.95 


0.95 


0.95 


0.3 


Linear 


0, 


,1 


CK = 


1 




0.98 


0.97 


0.96 


0.99 


0.97 


0.96 


0.8 


Linear 


0, 


,1 


a = 


1 




0.97 


0.95 


0.95 


0.99 


0.96 


0.95 



Note. — Table of xi foi' representative values of the simulated data. The 
output columns correspond to Galactic extinction law with Am fixed (1) and free 
(2), power law with Am fixed (3) and free (4), and linear law with Am fixed (5) 
and free (6). 



dent on the redshift of the lens, with more nearby lenses 
requiring stronger extinction to distinguish the different 
extinction laws. In general, we find that for A{V) > 0.3 
we can distinguish the different extinction laws, in par- 
ticular if we have constraints on Atti. The parameters 
of the power law fit tend to be more poorly constrained 
than those of the other fits. This is due to the added 
degree of freedom which can result in several sets of pa- 
rameters giving similar xt- We also see that if we can 
constrain the intrinsic magnitude ratio, the different ex- 
tinction laws become easier to distinguish. 

2.6.2. Extinction in both images 

We first simulate data where we apply the same ex- 
tinction law to both images (i.e., Ry and a are constant 
respectively) but at different strengths. Here we find 
that unless A'^{V) ~ A^{V) (i.e., in effect a very weak 
differential extinction), the fits converge to the starting 
parameters and give the shape of the 'real' extinction 



curves. 

To see what effect extinction in both images would 
have on a data set, we simulate data having different 
kinds of Galactic extinction laws, to find within which 
accuracy the input parameters of the more strongly ex- 
tinguished line of sight are found (see Fig. EJ. Again 
we find, that if the differential extinction between the 
images, E^{B -V) - E^{B -V), is low, then the fits 
are not well constrained and the resulting parameters 
represent the extinction of neither line of sight. How- 
ever, if one image is significantly more extinguished than 
the other then the parameters of the fit converge to the 
parameters of that line of sight as we expect (see dis- 
cussion in § 12.3(1 . It is not crucial that one image be 
non-extinguished, but it does need to have significantly 
lower extinction than the other line of sight. We find that 
for the fits to be within one sigma of the real parameters 
for the line of sight for the stronger extinction we need 

roughly E{B - V)weaker / E{B - V)stronger < 0.2 whcn 
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Fig. 6. — A sample plot for simulations with Galactic type ex- 
tinction in both images at two different redshifts (z = 0.0 and 0.8). 
The solid and the dotted lines are the input extinction laws and 
the dashed line is the effective differential extinction between the 
two images. We see that when the extinction along one line of sight 
is much stronger than the other, then the effective extinction law 
approaches the extinction of the more strongly extinguished image. 



the Rv values lie in the range of 2 — 4. This is consistent 
with our results from i)2.3l fin particular, see eq. (|14|l and 

2.6.3. The effects of noise 

When dealing with real data we can expect our data 
sets to be contaminated by various sources of noise. To 
see how this may affect our analysis we generate data 
sets with artificial random noise (see Fig. [7| for rep- 
resentative plots). We generate the noise as normally 
distributed random numbers with a mean of and stan- 
dard deviation of 0.05 magnitudes (which is an estimate 
of the lowest noise expected from deconvolved ground 
based data such as our data set). We still get quali- 
tatively the same results, i.e., we get the lowest for 
the Galactic type extinction fit to the Galactic extinc- 
tion law data (provided the extinction is strong enough) , 
and the fitted parameters agree with the input parame- 
ters, within the uncertainty. The values for the fits 
are increased and the uncertainty on the parameters for 
the noisy data is generally larger. If the noise becomes 
too large, i.e. of the order of the extinction effects, the 
fits become badly constrained. This emphasizes the need 
for a strong extinction signal to outweigh random noise 
effects when analyzing the real data sets. 

2.6.4. Achromatic microlensing 

We create several sets of data consistent with the 
Galactic extinction law with an effective microlensing pa- 
rameter 0.01 < s < 0.3 to test the effects of achromatic 
microlensing on our methods (see Fig. IHJ. We find that 
for very weak microlensing (s = 0.01) there is no no- 
ticeable effect on the results (given E{B — V)> 0.1), in 
particular not when we include noise (of 0.05 mag) in our 



Fig. 7. — A sample plot for simulations with noise showing three 
data sets with AiV) = 0.1,0.5 and 1.0. The diamonds and the 
dash-dotted line show the noiseless input extinction law (Galac- 
tic extinction law with Ry = 3.0), and the solid dots show the 
data points with random noise applied. The three fits shown are 
the Galactic extinction law (dashed line, eq. 0), power law (dot- 
ted line, eq. J7J) and linear law (solid line, |oJ). We see that for low 
extinction (A(V) < 0.1) the noise is of the order of the amount of 
the extinction signal, and therefore significantly affects the quality 
of the fits. 





0.8 
















0.6 
0.4 




s=0.30 










<: 

CD 


0.2 
0.0 
-0.2 
















-0.4 
















0.4 




5=0.20 






JiJ:-^- 




3-A (mag) 


0.2 
0.0 
-0.2 




I 




L 








-0.4 
















0.4 




s=0.10 










B-A (mag) 


0.2 
0.0 
-0.2 










i 






-0.4 
















0.4 




s=0.05 










B-A (mag) 


0.2 
0.0 
-0.2 
-0.4 
















0.0 


0.5 


1.0 


1.5 


2.0 2.5 


3 



.0 



Fig. 8. — A sample plot for simulations with an achromatic mi- 
crolensing signal s = 0.05,0.10,0.20,0.30 and noise. The asterisks 
are the noise free input data points, the boxes are the input data 
points with noise, the solid circles are the microlensed data points 
we fit and the diamonds are the microlensing corrected points from 
the fit. The dashed and dash-dotted lines are the best fitting Galac- 
tic extinction law with Am free and fixed through the microlensed 
data points and he solid line shows the best fitting Galactic ex- 
tinction law with Am free with a microlensing correction. The 
microlensing signal starts dominating the noise when the average 
distance from the asterisks to the boxes becomes smaller than the 
average distance from the boxes to the filled circles which happens 
around s = 0.05. 
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TABLE 2 
Table of Lenses 









iype 




Images 


"9 










doeidJ") 




(") 




Q2237+030 


0.04 


1.70 


Late'' 


-0.59 ±0.08 


{A,B,C,D) 


(0.92, 0.97, 0.76, 0.88) 


1, 2 


PG1115+080 


0.31 


1.72 


Early 


-0.33 ±0.02 


{A1,A2,B,C) 


(1.18, 1.12, 0.95, 1.37) 


1, 3 


B1422+231 


0.34 


3.62 


Early 


-0.50 ±0.13 


{A,B,C,D) 


(0.95,0.89,1.04,0.36) 


1, 4 


B1152+199 


0.44 


1.02 


Late'' 


-0.8±0.2<: 


(A,B) 


(1.14, 0.47) 


5, 6 


Q0142-100 


0.49 


2.72 


Early 


-0.29 ±0.02 


(A,B) 


(1.86, 0.38) 


1, 7 


B1030+071 


0.60 


1.54 


Early 


-0.35 ±0.06 


(A.B) 


(1.28, 0.58) 


1, 7 


RXJ0911+0551 


0.77 


2.80 


Early 


-0.17 ±0.04 


{A,B,C,D) 


(0.87, 0.97, 0.82, 2.24) 


8, 9 


HE0512-3329 


0.93 


1.57 


Late 


-1.0±0.3<= 


(A.B) 


(0.035, 0.66) 


10, 11 


MG0414+0534 


0.96 


2.64 


Early 


-0.11 ±0.08 


{A1,A2,B,C) 


(1.19, 1.17, 1.38, 0.96) 


1, 12 


MG2016+112 


1.01 


3.27 


Early 


-0.66 ±0.05 


(A.B) 


(2.48, 1.25) 


1, 11 



R eferences. — Clt IRusin et alJ <2003^. ('21 IRix et"an IT991), ( 31 IKristian et al.l tT9Ml. 

(4) 'Yce & EUingson " ^994"), (5) 'Mvcrs et al." fl99g). f 61 IRusin et all I2002al), (7) ILehar et alj 
(2000), (8) Rusinctal. I2002b), (9) Burud et al. 1199 83), (1 0) IGregg et alj 120001) . (11) CASTLES 
|http://www. cfa.harvard.edu/castles/ I, (12l iFalco et alj 119971) . 

Note. — The table lists various properties of the lensing systems known from the literature. The 
properties listed are the lens and quasar redshifts {zi and 2q), the type and the scale length ids) of the 
lensing galaxy, and the image names and their distance from the center of the lens galaxy {dg). The 
scale length, da, is the effective radius of a de Vaucouleurs profile fit from lHusin et al.l 12003) . The final 
column lists the references from which the values were obtained if they have not been previously quoted 
in the text. 
''The bulge is responsible for the lensing. 

''The spectra taken by IMvers et al.l 119991) shows O II emission line associated with the lens. We 
therefore type B1152+199 as a late type galax y. 

'^For the purposes of the plotting of Figure Iz!^ we adopt these scale lengths although they are not 
reported in the literature. 



data points as the noise dominates the effects of the mi- 
crolensing. For a stronger microlensing signal (0.05 < s) 
wc find that the effects are indeed noticeable but to be 
able to quantify them it is crucial to be able to constrain 
the intrinsic magnitude ratio difference. We find that a 
fit with a free intrinsic ratio can often give as good a fit 
(as measured by the xt) as the fit which allows for cor- 
rection due to the microlensing signal. This is because 
the effects of the microlensing can in part be mimicked 
by shifting the whole data set up and down along the 
magnitude difference axis by changing the intrinsic mag- 
nitude difference if the ratio of the continuum emission 
to the total emission is similar for the different bands 
(see eq. ((T7|l '). 

3. OBSERVATIONS 

The lens systems were chosen to fulfill the criteria that 
they have an image separation larger than 0'.'4 to ensure 
that the images of the quasar could be resolved, that they 
have declination 5 < 33° to be visible with the ESO Very 
Large Telescope (VLT) at Paranal observatory and that 
the lens and quasar redshifts be known in order to re- 
duce the number of unknowns when fitting the extinction 
curve. At the time of the application, this left us with 
ten systems, five doubly imaged quasars (doubles) and 
five quadruply imaged quasars (quads). The images are 
labeled according to the CfA-Arizona-Space-Telescope- 
LEns-Survey (CASTLES)^ notation (see Fig. EI). 

Multi waveband imaging observations of the 10 gravi- 
tational lens systems was obtained with the VLT. A list 
of the systems and their main properties known from the 
literature is given in Table [21 and a gallery of how they 

® http:/ /www. cfa.harvard.edu/castles/ 



appear in the VLT observations is shown in Figure El 
Optical observations (in the U, B, V, R, I and z band) 
was carried out with the FORSl instrument (which with 
the high resolution collimator has a pixel scale of O'.'l), 
and near infrared (NIR) observations (in the Jg, H and 
Ks band) were carried out with the ISAAC instrument 
(which has a pixel scale of 0"147). 

The data were collected in excellent seeing conditions 
(FWHM < 0'.'65, see FigureHUl) with mean seeing of 0'.'57 
for the full data set. Photometric conditions were not 
necessary since we are considering only relative photom- 
etry. An effort was made to carry out the different wave- 
band observations of each system as close in time as pos- 
sible to each other, to minimize the effects of time depen- 
dent intrinsic quasar variation and achromatic microlens- 
ing. For each system, the optical wavebands were ob- 
served on the same night in immediate succession, while 
the NIR observations were observed as close in time as 
scheduling allowed (the mean delay was 18 days, see Ta- 
ble O- The effects of time delayed intrinsic variations 
between the individual images is thus reduced to a pos- 
sible shift between the optical and NIR fluxes. Details of 
the observations are summarized in Table O 

3.1. Data Reduction 

The individual optical data frames were bias sub- 
tracted and flat fielded using the standard ESO pipeline, 
and cleaned for cosmic rays using Laplacian edge detec- 
tion ijvan Dokkuml20(n|) . In some of the frames the back- 
ground was found to have large-scale gradients over the 
field. To properly account for this, we ran all the frames 
through the SEXTRACTOR software ijBertin Sz Arnout j 
11996^) . with options set to save a full resolution interpo- 
lated background map, which was then subtracted from 
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TABLE 3 

Overview of observations 



Lens 


U 


B 


V 


R 


/ 


z 


Js 


H 




Delay 




(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(d) 


Q2237+030 


180 


60 


60 


40 


60 


40*= 




60 


60 


-13 


PG1115+080 


30 


15 


9 


9 


9 


ga 


180 


300 


180 


-20 


B1422+231 


180 


60 


40 


40 


40 


40*^ 




120 


120 


48 


B1152+199 


3000 


180 


80 


40 


60 


60*^ 


60 


60 


60 


47 


Q0142-100 


60 


60 


40 


40 


60 


60*^ 










B1030+071 


9300 


3600 


3000 


1400 


1600 


2000=" 


840 


1280 


2470 


-5 


RXJ0911+0551 


900 


30 


15 


15 


25 


15=" 


180 


216 


360 


5 


HE0512-3329 


30 


15 


9 


9 


9 


ga 


216 


216 


216 


-18 


MG0414+0534 




5400 


4200 


270 


60 


120=^ 


240 


240 


240 


-11 


MG2016+112 


2400 


1200 


720 


300 


400 


600'' 


480 


960 


1500 


-71 



Note. — Total exposure time (in seconds) and delay between optical and NIR observa- 
tions (in days), where a negative value denotes that the NIR were carried out before the 
optical. 

^The observing band was ^specia! • 
''The observing band was z (Gunn) 



Fig. 9. — Gallery of the 10 gravitational lensed quasars in the 
sample as they appear in the R-band VLT images. The size of the 
stamps is 6" X 6". 




0.0 0.5 1.0 1.5 
Wavelength (^un) 



2.0 



Fig. 10. — Mean FWHM measured in the VLT observations 
of the 10 systems, as a function of wavelength (full curve). The 
shaded area between the dashed curves indicates the RMS scatter 
around the mean. Approximately 80% of the observations were 
carried out in < 0.65" seeing. 



the science frames. 

The NIR data w ere reduced using a combination of the 
ECLIPSE software l|Devi llard 1997") and the IRAF "Ex- 
perimental Deep Infrared Mosaicing Software" xdimsum. 
The ECLIPSE software was used to remove the effects of 
electrical ghosts from science and calibration frames, and 
to construct flat fields and bad-pixel maps from a series 
of twilight sky flats. Sky subtraction and combination of 
the individual science frames was carried out with xdim- 
sum. 



3.2. Deconvolution and photometry 

All ph otometry was carried out using the MCS decon- 
volution ijMagain et al.ll99^ . This method uses a model 
PSF, measured directly from the data, to deconvolve the 
images, assuming that the gravitational lens systems can 
be decomposed into a number of point source (the quasar 
images) and a diffuse extended component (the lensing 
galaxy) . Positions and amplitudes of the quasar images 
are left free in the fit, and relative photometry can be 
derived from the best fitting amplitude. No functional 
form was assumed for the diffuse component, which is 
a purely numerical component. In the cases when there 
was more than one image of a given system available (in 
a given waveband), we performed simultaneous deconvo- 
lution of the individual images rather than deconvolving 
the combined image. 

Photometric errors were estimated by the MCS algo- 
rithm, and include p hoton noise and err o rs associated 
with deconvolution ijMagain et al.l Il998t iBurud et alJ 
[l998b). A fufl list of the results from the MCS decon- 
volution is given in Table 2| In addition to the MCS 
errors, we applied a 0.05 mag error on all the calculated 
magnitude differences to account for other sources of sys- 
tematic noise. 

We exclude data points from our sample where the 
deconvolution did not converge or where only one com- 
ponent was detected (7 data points). In addition, we ex- 
clude 4 data points from further analysis. These points 
are marked in Table ^ and an explanation of the exclu- 
sion of each point is given in a footnote. We also exclude 
from further analysis all the data taken for B1030-I-071 
as the B component is heavily cont aminated by the lens 
galaxy (separated by O'.'ll ± O'.'Ol ijXanthopoulos et all 
Il998i) ). 

4. RESULTS AND DISCUSSION 

We open this chapter by presenting in S 14. II the results 
of our extinction curve analysis for each of the 10 lensing 
systems. We present the detailed analysis of systems 
where at least one image pair has a two sigma detection 
of extinction for one of the three applied extinction laws. 
The systems are presented in order of increasing redshift. 



Extinction curves of lensing galaxies 



TABLE 4 



Lens 


image 


r 7 




V 


JX 


1 


z 


Js 


u 
ri 




Q2237+030 


A 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 




1.00 


1.00 






0.02 


0.02 


0.01 


0.01 


0.02 


0.02 




0.01 


0.01 




B 


0.36 


0.24 


0.30 


0.39 


0.30 


0.32 




0.33 


0.36 




a Q 


0.02 


0.02 


0.01 


0.01 


0.02 


0.02 




0.01 


0.01 




C 


0.26 


0.35 


0.40 


0.29 


0.37 


0.45 




0.47 


0.41 




(T<7 


0.02 


0.02 


0.01 


0.01 


0.02 


0.02 




0.01 


0.01 




D 


0.24 


0.21 


0.29 


0.27 


0.27 


0.35 




0.38 


0.35 




(7 ]j 


0.02 


0.02 


0.01 


0.01 


0.02 


0.02 




0.01 


0.01 


PG1115+080 


Al 


1.00 


1.000 


1.000 


1.000 


1.000 


1.00 


1.000 


1.000 


1.000 




(7 jW 


0.01 


0.006 


0.006 


0.005 


0.003 


0.01 


0.001 


0.002 


0.002 




A2 


0.63 


0.615 


0.676 


0.675 


0.710 


0.67 


0.747 


0.710 


0.771 






0.02 


0.009 


0.007 


0.008 


0.004 


0.02 


0.002 


0.003 


0.002 




B 


0.26 


0.246 


0.253 


0.257 


0.262 


0.27 


0.268 


0.270 


0.266 






0.02 


0.009 


0.008 


0.008 


0.004 


0.02 


0.002 


0.003 


0.003 




C 


0.18 


0.17 


0.17 


0.17 


0.170 


0.17 


0.177 


0.178 


0.176 






0.02 


0.01 


0.01 


0.01 


0.006 


0.01 


0.003 


0.005 


0.005 


B1422+231 


A 


0.105 


0.15 


0.16 


0.15 


0.18 


0.16 






0.31=^ 






0.008 


0.02 


0.01 


0.01 


0.02 


0.01 






0.02=^ 




B 


0.09 


0.14 


0.16 


0.15 


0.17 


0.16 






0.35=^ 




ag 


0.01 


0.02 


0.02 


0.02 


0.02 


0.01 






0.02=^ 




C 


0.08 


0.10 


0.10 


0.09 


0.11 


0.10 






0.18=^ 




""C 


0.02 


0.03 


0.04 


0.03 


0.04 


0.01 






0.02=^ 




D 


0.01 


0.004 


0.005 


0.005 


0.007 


0.006 






0.02=^ 






0.01 


0.018 


0.016 


0.014 


0.019 


0.013 






0.02=^ 


B1152+199 


A 




1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 








0.007 


0.002 


0.002 


0.003 


0.005 


0.004 


0.002 


0.004 




B 




0.0024 


0.010 


0.018 


0.044 


0.074 


0.172 


0.293 


0.289 








0.0006 


0.001 


0.001 


0.002 


0.003 


0.002 


0.004 


0.003 


C4U14Z— lUU 


A 

A 


1.000 


1.00 


1.0000 


1.0000 


1.0000 


1.0000 










(^A 


U.UUO 


U.Ul 


U.UUUO 


n 0009 


u.uuuz 


U.UUUO 










B 


0.13 


0.12 


0.134 


0.138 


0.152 


0.159 










(7 b 


0.03 


0.01 


0.001 


0.001 


0.001 


0.002 








B1030+071^ 


A 




1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 




^ A 




0.01 


0.01 


0.01 


0.01 


0.02 


0.02 


0.03 


0.02 




B 




0.16 


0.23 


0.30 


0.41 


0.64 


0.378 


0.29 


0.28 








0.01 


0.01 


0.01 


0.01 


0.02 


0.008 


0.01 


0.01 






1.000 


1.00^ 


1.00 


1.00 


1.00 




1.000 


1.000 


1.000 






n nnc 


n nob 


U.Ui 


U.Ui 


n no 
U.Uz 




U.UUo 


U.UUo 


U.UUO 




B 


0.62 


0.94'' 


0.73 


0.72 


0.74 




0.923 


0.919 


0.970 




o-s 


0.01 


0.02'' 


0.02 


0.02 


0.02 




0.009 


0.008 


0.006 




C 


0.33 


0.51'' 


0.37 


0.40 


0.41 




0.49 


0.51 


0.50 






0.02 


0.03*^ 


0.04 


0.03 


0.04 




0.02 


0.01 


0.01 




D 


0.37 


0.40*^ 


0.38 


0.37 


0.37 




0.39 


0.44 


0.40 






0.01 


0.02'' 


0.02 


0.01 


0.02 




0.02 


0.02 


0.01 


HE0512-3329 


A 


1.000 


1.000 


1.000 


1.0000 


1.000 




1.0000 


1.0000 


1.0000 






0.008 


0.002 


0.002 


0.0009 


0.001 




0.0005 


0.0004 


0.0005 




B 


1.141 


1.021 


0.887 


0.738 


0.651 




0.5813 


0.5667 


0.5593 




o-s 


0.007 


0.002 


0.002 


0.001 


0.002 




0.0008 


0.0007 


0.0009 


MG0414+0534 


Al 






1.000 


1.000 


1.00 


l.OO'^ 


1.000 


1.000 


1.000 










0.005 


0.008 


0.01 


0.02=^ 


0.004 


0.003 


0.001 




A2 






0.41 


0.34 


0.39 


0.89=^ 


0.567 


0.748 


0.780 




<^A2 






0.01 


0.02 


0.03 


0.02== 


0.008 


0.004 


0.003 




B 






0.818 


0.592 


0.50 


0.78^^ 


0.416 


0.394 


0.385 










0.005 


0.009 


0.01 


0.02=: 


0.006 


0.004 


0.006 




C 






0.396 


0.29 


0.24 


0.34^= 


0.18 


0.179 


0.175 










0.008 


0.01 


0.02 


0.03^= 


0.01 


0.006 


0.005 


MG2016+112 


A 




1.00 


1.00 


1.00 


1.00 


l.OO'' 


1.0 


1.00 










0.01 


0.01 


0.01 


0.01 


0.02'' 


0.1 


0.03 






B 




0.60 


0.68 


0.81 


0.87 


1.04'' 


1.08 


0.94 










0.02 


0.01 


0.01 


0.01 


0.02^ 


0.09 


0.03 





Note. — Data table with tiie results from the deconvolution. Missing data are either due to lack of 
observations (see overview in Table HI or failure of the deconvolution to converge. 

tThe photometry of the B component is believed to be contaminated by the lensing galaxy. 
^The images are not well resolved. This system shows very weak extinction, and no detailed extinction 
curve analysis is performed, so the exclusion of this point should not affect our results. 
''Seeing too high to separate components (0''70). 
"^Very faint detection. 
''Very faint detection. 
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We then move on to discussing in i) 14 .21 the overaU resuhs 
of our analysis, as weU as statistical properties of the full 
sample. 

4.1. The individual systems 
4.1.1. Q 2237+030 

Q2237+030 was discovered by iHuchra et all l|1985fl 

and consists of a quasar at redshift z = 1.70 and a 
spiral lensing galaxy at z = 0.04 making it the nearest 
known lensing galaxy to date. The system was later re- 
solved into four images formin g an Einstein cross with the 
lensing galaxy in the middle (lYeeHlQSSt iSchneider et al.l 
^88). Schneider et al. (198^ modeled the system in de- 
tail and found a predicted time delay of order of one day 
between the images and amplifications of 4.6, 4.5, 3.8 and 
3.6 fo r images A, B, C and D, respectively. iFalco et aP 
l)1996j) studied the system with the VLA at radio wave- 
lengths and obtained flux density ratios of 1.00, 1.08, 0.55 
and 0.77 for images A, B, C, D, respectively. Q2237+030 
has previously been noted in the literature as having high 
variabilit y which is uncorrelated between the fou r images 
(see e.g.. llrwin eraLlll989tlCorrigan et alJll99lD . 

This system is difficult to interpret as it shows a lot 
of scatter in the points which cannot be explained by 
extinction alone nor microlensing (see Figures 1111 and 
112(1 . None of the extinction laws we apply give a good 
fit to the data for any pair of images (as the redshift of 
this system is very low, z = 0.04, we do not expect to 
see much difference between the quality of the fits for 
the different extinction laws, see S I2.6.1|) . All the image 
pairs, except C— A and D— A, yield A{V) consistent with 
zero. The data points and fits for C— A and D— A can be 
seen in Figures ^2 and El and the parameters of the fits 
in Tables El and ini 

We use the radio measurements of IFalco et al.l 1)19961) 
to fix Am and in particular for the D— A image pair this 
changes the results significantly. As the flux density ra- 
tios in the radio agree with the model predictions for the 
D— A image pair it is interesting to note that none of 
the extinction laws we apply give good fits to the data 
unless we allow for corrections due to achromatic mi- 
crolensing (see Figures ^2 and ^1 and Tables [SI and 0) . 
This might suggest that either D is dcmagnificd by a 
strong microlensing signal, or that image A is magnified 
(or both) and the residual 'extinction curve' which we 
are fitting may be effects of chromatic microlensing. The 
same effect, but not as strong, is seen in the C— A image, 
again suggesting a slight demagnification of C, a magni- 
fication of A or a combination of the two. In previous 
microlensing studies, component D has not been seen to 
have as strong a microlensing signal as A and C have 
lllrwin et al.lll989l: IAlcaldell200l IGil-Merino eTalWm)^ 
so therefore it is perhaps more likely that we are seeing 
the magnification of image A. 

Another explanation for the shift in Am might be in- 
trinsic variations of the quasar components as discussed 
in '<i2.1\ Such variability could introduce an overall shift 
of the data, resulting in inaccurate estimates of Am. 

By inspecting the l^-band lightcurves from iKochanekl 
(|2"004) we see that component A is indeed in a bright 
phase at the time of our observations (Julian date of 
around 2451670). Component D is fairly stable but com- 
ponent C is getting dimmer climbing down from a peak 




2 4 6 

V (urn-') 

Fig. 11.— Q2237+030, C-A: The upper panel shows the data 
points and the best fit extinction curves. The lower panel shows 
the original data points and their shift due to a microlensing signal, 
along with the best fit through the shifted points. The parameters 
of the fits are given in Tablel^ Annotation: Filled circles are the 
original data points with error bars. The curves correspond to the 
Galactic extinction law (e g. JE\ with Am free (dash-dot) and fixed 
(solid), the power law (eq. 171 with Am free (dash-dot-dot-dot) and 
fixed (dotted) and the linear law (eq. |SJ with Am free (long dash) 
and fixed (short dash). Shifted data points due to a microlensing 
signal are plotted in open boxes (Galactic extinction law), triangles 
(power law) and diamonds (linear law). 




cot . . . : . . . : . . . J 
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Fig. 12.— Q2237+030, D-A: The upper panel shows the data 
points and the best fit extinction curves. The lower panel shows 
the original data points and their shift due to a microlensing signal, 
along with the best fit through the shifted points. The para met ers 
of the fits are given in Table l6l See the caption of Figure ITTl for 
annotation overview. 



in its brightness and is still fairly bright, perhaps making 
the C— A shift in Am less prominent than the D— A shift. 

4.1.2. PG1115+080 



Extinction curves of lensing galaxies 
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TABLE 5 

Extinction curve fit results for Q2237+030: C— A 



Extinction 


A{V) 


Rv 


a 


Am 


s 




Galactic 


0.53 ±0.03 


> 7 




0.65 




3.6 ±0.4 


Galactic 


0.29 ±0.06 


2.9 ± 1.4 




0.85 ±0.04 




3.4 ±0.4 


Galactic 


0.27 ±0.07 


2.6 ± 1.5 




0.65 


-0.23 ±0.05 


3.4 ±0.4 


Power law 


0.49 ±0.02 




0.6 ±0.1 


0.65 




3.2 ±0.4 


Power law 


0.3 ±0.5 




1±4 


0.8 ±0.5 




3.4 ±0.4 


Power law 


0.3 ±0.2 




1±3 


0.65 


-0.2 ±0.1 


3.4 ±0.4 


Linear law 


0.49 ±0.02 




1.0 


0.65 




3.2 ±0.4 


Linear law 


0.35 ±0.05 




1.0 


0.78 ±0.04 




3.1 ±0.4 


Linear law 


0.34 ±0.05 




1.0 


0.65 


-0.15 ±0.05 


3.1 ±0.4 



Note. — The extinction of image C compared to image A. Numbers quoted in italics 
were fixed in the fitting procedure. 



TABLE 6 

Extinction curve fit results for Q2237+030: D— A 



Extinction 


A{V) 


Rv 


a 


Am 


s 


xl 


Galactic 


1.23 ±0.03 


> 7 




0.28 




7.1 ±0.4 


Galactic 


0.35 ±0.06 


3.1 ± 1.5 




1.05 ±0.05 




1.9 ±0.4 


Galactic 


0.28 ±0.07 


2.1 ± 1.2 




0.28 


-0.85 ±0.05 


1.9 ±0.4 


Power law 


1.11 ±0.03 




0.27 ±0.04 


0.28 




1.9 ±0.4 


Power law 


0.4± 1.1 




0.9 ±0.4 


0.9 ± 1.1 




2.0 ±0.4 


Power law 


0.3 ±0.4 




1.3 ±0.6 


0.28 


-0.8 ± 0.4 


2.0 ±0.4 


Linear law 


1.18 ±0.03 




1.0 


0.28 




5.9 ±0.4 


Linear law 


0.41 ±0.06 




1.0 


0.97 ±0.05 




1.9 ±0.4 


Linear law 


0.39 ±0.06 




1.0 


0.28 


-0.75 ±0.05 


1.9 ±0.4 



Note. — The extinction of image D compared to image A. Numbers quoted in italics 
were fixed in the fitting procedure. 



PG1115+080 is a multiply imaged system discovered 
by iWevnmann et a l. (1980) as a triply imaged system 
with the quasar at redshift z = 1.72. The A compo- 
nent wa s later resolved int o two separate images, Al and 
A2, by iHeee et all 1)19811) maki ng the system a quad . 
The lensing galaxy was loca ted bvlChristian et alJ l|1987f) 
and is an early type galaxy ijRusin et al.ll2002atl . Its red- 
shift and that of three neighboring galaxies were deter- 
mined to be at z = 0.31 by Kundic ct al. (1997a). The 
time delays between th e components were determined by 
ISchechter et"all l)1997|) . 

The system shows very weak differential extinction for 
all the images (all but A2— Al have A{V) equal to 
within two sigma, see Tabic ITT)). The data points and fits 
for A2— Al are shown on Figure ITSl and the parameters 
of the fits in Table [7| The low extinction sig nal is in 
agreement with the results of iFalco et all l)1999|) . 

4.1.3. BU22+231 

B1422+2 31 is a quadruply irn aged system first dis- 
covered by IPatnaik et all l)1992|) in th e JVAS survey 
and co nfirmed to be a lensing system bv lLawrence et all 
()1992j) . The l ensing system consists of an early type 
main galaxv llYee Ellingsonl 11994 at z = 0.34 
llK nndic et al.ll9971: . Tonrvl 19981) and five n earbv galax- 
ies ijRemv et al.ll993tlBechtold fc Yeell995|) . The quasar 
is at a redshift of z = 3.62 and the m aximum image 
separation is 1.3" ('Patnai k et "all 119921) . The images 
show intrinsic variability which has been used to de- 
termine the time delay by studying radio light curves 




2 4 



Fig. 13.— PG1115+080, A2-A1: Best fit extinction curves for 
A2— Al. The param eter s of the fits can be seen in Table 171 See 
the caption of Figure [TTl for annotation overview. 



ijPatnaik k. Narasimhall2001|) . 

We only use three of the four images in our analysis 
as D was too faint to give usable results (all the visi- 
ble bands gave zero detection). As for the other com- 
ponents they show very weak differential extinction and 
give very weak constraints on the differential extinction 
curves. All the fits have AiV) consistent with zero (see 
Table ITTj) . We also get AiV) consistent with zero when 
we fix Am (where the values for the A m are taken to 
be the average between those deduced by IPatnaik et alJ 
lfT99^ in the 5 GHz and 8 GHz bands). The low dif- 
ferential extinctio n between the ima ges is in agreement 
with the results of IFalco et alJ l)1999D . 
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TABLE 7 

Extinction curve fit results for PG1115+080 



Extinction 


A{V) 


Rv 


a Am 




Galactic 


0.12 ±0.06 


3.3 ± 1.9 


0.29 ±0.05 


1.2 ±0.3 


Power law 


0.3 ±2.0 




0.5 ±1.5 0.1 ±2.0 


1.2 ±0.3 


Linear law 


0.13 ±0.03 




1.0 0.28 ±0.04 


1.2 ±0.3 



Note. — The extinction of image A2 compared to reference image Al 
Numbers quoted in italics were fixed in the fitting procedure. 



4.1.4. B1152+199 

B1152+199 is a doubly imaged system first discovered 
by M yers et al... (1999) in the CLASS survey with a back- 
ground quasar at z = 1.02, a lensing galaxy aX z — 0.44 
and image separation of 1'.'56. It was observed in radio 
wavelengths (at fre quencies 1.4,5,8.4 and 15 GHz) by 
iRusin et aLll)2002aj) . The extinction curve has previously 
been studied and fitted by a Galact ic extinction law with 
1.3 <Rv < 2.QaTulE{B-V) ~ 1 (|Toft et al.l200nD sug- 
gesting that it is a heavily extinguished system. 

B1152+199 shows a very strong extinction signal as 
can be seen in Figure^land Table|Sl It has the strongest 
extinction signal of all ten systems with AiV) = 2.43 ± 
0.09, 2.7 ±0.1, 3.57 ±0.07 at xl = 2.0,2.8,3.4 for the 
Galactic extinction law, power law and li near law re- 
spectiv ely. Using the radio measurements of lR.usin et al.l 
(j2002a(l to fix Am we similarly get A{V) = 2.03 ± 
0.03, 2.01 ± 0.03, 1.94 ± 0.03 at xl = 2.7, 4.0, 10.0. We 
also analyze the data with respect to a possible achro- 
matic microlensing signal, keeping Am fixed. This yields 
a non-zero microlensing correction for the Galactic ex- 
tinction law and power law of s = 0.32 ± 0.07, 0.6 ± 0.1 
and A{V) = 2.41 ± 0.09, 2.6 ± 0.1 at xl = 2.1,3.2 re- 
spectively. The best fit for the linear law lies outside the 
validity of the method with s > 1 (see i) l2.4|l which would 
correspond to a microlensing signal of > 1 mag. 

It is clear that in all cases the Galactic extinc- 
tion law provides the best fit to the data suggesting 
Galactic type dust although the best fit Rv values are 
lower than those commonly seen in the Milky Way. It 
is possible that the measured Rv value is being lowered 
by a non-zero extinction in the A image provided it has 
a higher value of Rv (see discussion in S 12.3(1 . However, 
given the very strong extinction signal this would require 
very strong extinction along both lines of sight in addi- 
tion to a strong differential signal. This is unlikely given 
the fact that component A is at more than twice the 
distance from the center of the lensing galaxy than com- 
ponent B with A at 1'.'14 and B at 0'.'47 from the center 
IIK.usm e i~ani20n2all . Measurements in the X-ray further 
suggest that the A component is non-extinguished (K. 
Pedersen et al., 2006, in preparation). 

4.1.5. QOI42-IOO 

Q0142— 100 is a dou bly imaged system first discovered 
bv lSurdei et al.l ^^1987^ and also known in the literature 
as UM 673. The quasar is at a redshift of z = 2.72 
(^lac Alpine & Fcldman 198 ^ and the le nsing galaxy, 
which is of early t ype (Rusin ct al. 2002a), is at a red- 
shift of z = 0.49 (|Su"rdei et al. 1987) . Wisotzki et all 
l|2004D studied this system using spectrophotometric ob- 
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Fig. 14. — B1152+199: The upper panel shows the data points 
and the best fit extinction curves to them. The lower panel shows 
the original data points and their shift due to a microlensing signal. 
The p aram eters of the fits are given in Table See the caption of 
Figure [TTl for annotation overview. 



servations and found signs of differential extinction but 
no microlensing. 

The data points and our fits for Q0142— 100 can be seen 
in Figure El and the parameters of the fits in Table El 
All the fits give similar xl but the parameters, in partic- 
ular for the power law, are poorly constrained due to the 
lack of data points (we did not get any measurements 
in the infrared for this system). The extinction is high 
for an ea rly type galaxy an d is not consistent with that 
found bv lFalco et alJ l|1999D who found negligible extinc- 
tion. We suspect that the data are being contaminated 
by the lens galaxy as the B component is located near 
the galaxy center (at 0'.'38) and the seeing was not opti- 
mal for this system (the mean seeing was 0'.'87 compared 
to 0'.'57 for the full data set). As there are no published 
radio measurements available we do not have constraints 
on Am to analyze the system with respect to a possible 
microlensing signal. 

4.1.6. B1030+071 

B1030±071 is a doub l y ima ged system first discovered 
bv lXanthopoulos et alJ (l998D in the JVAS survey. They 
monitored the system in radio wavelengths, finding that 
the flux density ratios between image A and B range 
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TABLE 8 

Extinction curve fit results for B1152+199 



Extinction 


A{V) 


Ry 




a 




Am 


s 




Galactic 


2.03 ±0.03 


1.61 ±0.05 








1.18 




2.7 ±0.4 


Galactic 


2.43 ± 0.09 


2.1 ± 0.1 








0.85 ± 0.07 




2.0 ±0.4 


Galactic 


2.41 ±0.09 


2.0 ±0.1 








1.18 


0.32 ±0.07 


2.1 ±0.4 


Power law 


2.01 ±0.03 




1.98 


± 


0.04 


1.18 




4.0 ±0.4 


Power law 


2.7 ±0.1 




1.45 


± 


0.08 


0.6 ±0.1 




2.8 ±0.4 


Power law 


2.6 ± 0.1 




1.52 


± 


0.07 


1.18 


0.6 ±0.1 


3.2 ±0.4 


Linear law 


1.94 ±0.03 








1.0 


1.18 




10.0 ±0.4 


Linear law 


3.57 ±0.07 








1.0 


-0.23 ± 0.06 




3.4 ±0.4 



Note. — The extinction of image B compared to reference image A. Numbers quoted in 
italics were fixed in the fitting procedure. 



TABLE 9 

Extinction curve fit results for Q0142— 100 



Extinction 


A{V) 


Rv 


a 


Am 




Galactic 


0.40 ±0.03 


4.8 ±0.7 




1.64 ±0.04 


1.2 ±0.5 


Power law 


4.1 ± 3.8 




0.1 ±0.8 


-2.1 ± 3.8 


1.1 ±0.5 


Linear law 


0.27 ±0.08 




1.0 


1.8 ±0.1 


1.1 ±0.4 



Note. — The extinction of image B compared to reference image A. 
Numbers quoted in italics were fixed in the fitting procedure. 



iBurud et al.l l)1998af) who resolved the system into four 
images and found that large external shear, possibly due 
to a clus ter, was re quired to explain the image config- 
uration. iKneib et~a l. (2000) confirmed that the lensing 
galaxy belongs to a cluster dX z = 0.769. Observed red- 
dening in at least two (images B and C) of the four images 
suggest differential extinction by the early typ e lensing 
galaxy (Burud et al. 1998a). iHiorth et all l|2iD02il mea- 
sured the time delay of the system between images A,B,C 
on the one hand and D on the other and found the time 
delay to be 146 ± 8 days {2a). 

We find relatively strong extinction in images B and 
C compared to images A and D. Image D also shows 
some extinction when compared to A but the effect is 
consistent with zero within two sigmas. We analyze the 
extinction curves of B and C compared to A and D. The 
data points and the fits can be seen in Figures IT?^ and [T7I 
and the parameters of the fits can be seen in Tables [TUl 
and HI] 

Assuming image A is completely unextinguished we 
can estimate the lower limit of the relative extinction of 
D compared to B and C, Ed{B - V)/Eb,c{B - V). For 
both B and C, we find that this ratio is around 1/3 so we 
expect the extinction curve properties to be affected by 
both lines of sight (see discussion in i) l2.3|) . That is, we do 
not expect the extinction curve we get from comparing 
B and C to D to represent the extinction curve along 
either line of sight unless their extinction properties are 
identical. We see that the Ry value of both B and C 
are lower when compared to image D than those we get 
from comparing them to image A suggesting that the 
extinction properties are indeed different (with image D 
having a higher Ry value). We note however that the 
values of Ry do agree within one sigma for both the 
differential extinction curves for both B and C. 
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Fig. 15.— Q0142-100: The plot shows the data points and the 
best fit to them. The par ame ters of the fits can be seen in Tabled 
See the caption of Figure ITTl for annotation overview. 



from 12.0 to 18.8 and seem to vary with both time and 
frequency. The redshift of the background source was 
determined to be at z = 1.54 and the redshift of the lens- 
ing object to be a t z = 0.60 ijFassnacht fc Co hen 1998). 
iFalco et al.l ()1999f l determined a differential extinction of 
E(B -Y) = 0.02 ± 0.04 assuming a fixed Ry = 3.1 
Galactic extinction law. 

We were unable to perform an extinction analysis on 
this system as the deconvolution did not succeed in sepa- 
rating the B compon ent from the main lens gala xy (sepa- 
rated by O'.'lliO'.'Ol ijXanthoDoulos et al.ll998D ) making 
the photometric values unreliable. For a further study of 
the extinction of this system higher resolution images 
would be required. 

4.1.7. RXJ0911+0551 

RXJ0911+0 551 is a multiply imaged system first dis- 
covered bv .Bade et al.l (1997 ) in the ROSAT All-Sky Sur- 
vey with the quasar at z = 2.80. It was later studied by 
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TABLE 10 

Extinction curve fit results for RXJ0911+0551: B 



Image pair 


Extinction 


A{V) 


Rv 


a 


Am 




B-A 


Galactic 


0.33 ±0.06 


4.9 ±0.6 




-0.10 ±0.05 


1.3 ±0.4 


B-D 


Galactic 


0.20 ±0.08 


4.1 ± 1.1 




-1.00 ±0.06 


1.6 ±0.5 


B-A 


Power law 


0.9 ± 3.7 




0.3 ±0.3 


-0.7± 3.7 


1.3 ±0.4 


B-D 


Power law 


0.2 ± 1.4 




0.8 ±0.6 


-1.1 ± 1.4 


1.6 ±0.5 


B-A 


Linear law 


0.23 ±0.03 




1.0 


-0.04 ±0.04 


1.5 ±0.4 


B-D 


Linear law 


0.17 ±0.03 




1.0 


-1.00 ±0.04 


1.5 ± 0.4 



Note. — The extinction of image B compared to reference images A and D. Numbers 
quoted in italics were fixed in the fitting procedure. 



TABLE 11 

Extinction curve fit results for RXJ0911+0551: C 



Image pair 


Extinction 


A(V) 


Ry 


a 


Am 


xl 


C-A 


Galactic 


0.20 ±0.10 


3.3 ± 1.5 




0.66 ±0.07 


1.2 ±0.4 


C-D 


Galactic 


0.09 ± 0.08 


2.2 ± 1.5 




0.26 ±0.06 


1.3 ±0.4 


C-A 


Power law 


0.3 ± 2.1 




0.9 ±0.5 


0.6 ±2.1 


1.2 ±0.4 


C-D 


Power law 


0.1 ± 1.8 




1.3±1.5 


03 ± 1.8 


1.2 ±0.4 


C-A 


Linear law 


0.23 ± 0.04 




1.0 


0.62 ±0.04 


1.1 ±0.4 


C-D 


Linear law 


0.17 ±0.04 




1.0 


0.31 ±0.05 


1.2 ±0.4 



Note. — The extinction of image C compared to reference images A and D. Numbers 
quoted in italics were fixed in the fitting procedure. 
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Fig. 16. — RXJ0911+0551, B: The upper panel shows the data 
points and the fits for the extinction of image B compared to im- 
age A. The lower panel shows the corresponding plot for image B 
comp ared to image D. The parameters of the ilts can be seen in 
Table ITol See the caption of Figure ITTl for annotation overview. 



Fig. 17.— RXJ0911+0551, C: The upper panel shows the data 
points and the fits for the extinction of image C compared to im- 
age A. The lower panel shows the corresponding plot for image C 
comp ared to image D. The parameters of the ilts can be seen in 
Table ITTl See the caption of Figure ITTl for annotation overview. 



4.1.8. HE0512~3329 

HE0512— 33 29 is a doubly imaged system first discov- 
ered by Gregg et al.l l|200Q.l with an image separation 
of 0'.'644 and quasar redshift of z = 1.565. They es- 
timated a redshift of z = 0.9319 for the lensing ob- 
ject and found that the lens is most likely a spiral 
galaxy. In addition, they estimate the differential red- 
dening assuming negligible microlensing and a standard 
Galactic extinction law with Ry = 3.1. This yields 



A{V) = 0.34 with the A image being redder than the 
B image. iWucknitz et alJ (|2003') worked further on dis- 
entangling microlensing and differential extinction, and 
estimated A{V) — 0.07 with A being the extinguished 
image. This fit results in an effective Ry = —2.0 which 
can be achieved if the two lines of sight have different 
Rv ■ 

In the case of HE0512— 3329, it is the brighter image, 
A, which shows extinction with respect to the B image. 
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TABLE 12 

Extinction curve fit results for HE0512— 3329 



Extinction 


A{V) 


Rv OL 


Am 


xi 


Galactic 


0.14 ±0.04 


1.7 ±0.4 


-0.67 ±0.03 


2.1 ± 0.4 


Power law 


0.23 ±0.09 


1.3 ±0.3 


-0.76 ± 0.09 


1.4 ±0.4 


Linear law 


0.35 ±0.02 


1.0 


-0.86 ± 0.04 


1.4 ±0.4 



Note. — The extinction of image A compared to image B. Numbers quoted 
in italics were fixed in the fitting procedure. 
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Fig. 18.— HE0512-3329: The plot shows the data points and 
the b est fit to them. The parame ters of the fits can be seen in 
Table IT^ See the caption of Figure [TTl for annotation overview. 



The system is interesting as one of the redshifted data 
points falls in the range where the 2175 A bump in the 
Galactic extinction law should lie (see Figure [TH|l . There 
is however no sign of a bump at A = 2175 A and both 
the power law and the linear extinction law give a much 
better fit (see Table ITS] for the parameters of the fits). 
We redo the fits with no constraints on the Rv values to 
see if our data could be fit by a negative Rv value but 
this does not change the result of Rv = 1-7 ± 0.4. As 
there is no radio data available we do not constrain the 
intrinsic ratio in the fits and we can not constrain the 
microlensing signal. 

Our results are not in agreement with those of 
IWuck nitz ct al. (20Q3) who found that fits with the 2175 
A bump better reproduced their data than those with- 
out, although the result was not highly significant. In 
addition, they found that it is crucial to take microlens- 
ing into account when analyzing the extinction curve, 
which might explain the discrepancy. However, the de- 
tected microlensing signal is only important at wave- 
lengths lower than those we probe, with a small possible 
effect in the B- and F-bands. Therefo re, a microlens- 
ing signal consistent with the results of 'Wuc knitz et ahl 
(p003) should not affect our results significantly. In addi- 
tion, we note that for their best fitting Rv their fit curves 
downwards for < 1 ^m~^ which is not consistent 
with our measurement in the K-hand (see Figure [TH|l . 

4.1.9. MGO414+O534 

MG0414-I -0534 is a quadruply imaged system first dis- 
covered by iHewitt et all l)1992() with image separation 
of up to 2". The quasar, at redshift oi z = 2.64, 
shows evidence of b eing heavily reddene d by dust in 
the lensing galaxy ijLawrence et all Il995f) . The lens, 
which has early type spectrum, is at redshift z = 



0.9584 llTonrv fc Kocha^lTool and was modeled by 
iFalco et al.l 1)19971) who found the brightness profile to 
be well represented by a de Vaucouleurs profile which is 
characteristic of an elliptical galaxy. Falco ct al. (199j|) 
studied the extinction curve of this system and fitted it 
to a Galactic extinction law giving a best fit of Rv — 
1.5, assuming that all lines of sight have the same Rv ■ 
lAngonin-Willaime et al.l l|199 9^ studied the origin of the 
extinction and found, that while the differential extinc- 
tion is likely due to the lensing galaxy, then there is also 
evidence fo r significant red dening which is intrinsic to 
the source. iKatz fc HewittI (0,993) did an extensive radio 
survey of the system and found that there was no sign 
of variability in the radi o flux ratios between their mea- 
surements and those of IHewitt et al.l l)1992D except for 
the C/B image ratio. 

For MG0414+0534, components Al and A2 show ex- 
tinction when compared to images B or C with A2 being 
the more strongly extinguished image (see TablelT7| . For 
the Al image we find different effective extinction laws 
depending on whether we compare with image B or C 
(see Table El and Figure IT^ . In all cases the power 
law gives the best fit and the linear law t he worst (we 
use the radio measurements of IKatz fc Hew itt ( 1993) to 
fix Am). For the Galactic extinction law the Rv val- 
ues do not agree suggesting that perhaps the extinction 
of images B and C is affecting the differential extinction 
curve. We also note though, that we would expect A{V) 
for Al— B to be around 0.3, to be consistent with the 
other values in Table El but the best fitting values give 
a lower value. We therefore perform another fit where we 
fix A{V) = 0.3 in the fits for the Al— B pair and this gives 
Rv — 2.1 ± 0.2; 1.8 ± 0.1 for Am fixed and free respec- 
tively, which are marginally consistent with the results 
compared to the C image. However, the xt = 2.6; 2.0 of 
these fits are significantly worse than those of the orig- 
inal fits. We do not see any evidence for microlensing 
except in the case of SMC-like linear extinction which 
still results in a worse fit than the other two extinction 
laws. 

For image A2 the Galactic extinction law gives the 
best fit when Am is kept fixed but otherwise the dif- 
ferent extinction laws give similar results (see Table ITTI 
and Figure EH). The parameters of the Galactic extinc- 
tion law are consistent when compared with images B 
and C suggesting that either A2 dominates the extinction 
signal or that B and C have similar extinction properties. 
There is no evidence for microlensing except in the case 
of the linear extinction law. We note that the absolute 
extinction of image A2, which must be greater or equal 
to the differential extinction in Table 1141 i s verv h igh 
given that the lens is an early type galaxy l)Goudfrooil 
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Fig. 19.— MG0414+0534, Al: The upper panel shows the data 
points and the fits for the extinction of image Al compared to 
image B. The lower panel shows the corresponding plot for image 
Al compared to image C. The param eter s of the fits can be seen 
in Table [T^ See the caption of Figure [TTI for annotation overview. 
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Fig. 20.— MG0414+0534, A2: The upper panel shows the data 
points and the fits for the extinction of image A2 compared to 
image B. The lower panel shows the corresponding plot for image 
A2 compared to image C. The parameters of the fits can be seen 
in Table IT^ See the caption of Figure [TTI for annotation overview. 



l)1994|) find A{V) < 0.35 for their sample of early type 
galaxies) . 

As the extinction of Al is significant compared to A2 
we expect the extinction properties of both lines of sight 
to affect the A2— Al extinction curve (see S I2.3|l . The 
fit of A2— Al for the Galactic extinction curve gives us 
Rv = 3.5 ± 0.4, 3.8 ± 0.7 at A(y) = 0.53 ± 0.03, 0.61 ± 
0.11 when Am is kept fixed or free respectively. If we 
assume that images B and C have zero extinction we 
can calculate the effective Ry we expect to get from eq. 



(O. The results can be seen in Table [Tsl and are in good 
agreement with the results of the fits. 

The extinction of MG0414+0534 is high for an early 
type galaxy. We can not exclude the possibility that the 
extinction may be due to an unknown foreground object 
and not the lensing galaxy itself. Finally we note that our 
es timates of the diffe rential extinction agree with those 
of iFalco et all l)l999j) which were obtained by assuming 
standard Galactic extinction with Ry =3.1. 

4.1.10. MG2016+112 

MG 2016+112 was discovered by iLawrence et alJ 
l)1984|) and ha s a giant ellip tical lensing galaxy at red- 
shift z = 1.01 l|Schneider et al. 1985, 1986). The system 
consists of two images, A and B, of the quasar at red- 
shift z — 3.273 and an additional image C which may 
be a third image of the quasar with an additional signal 
from another ga laxy and has been challenging to model 
(ILawrence eral. 1984. 1993: Nair fc Garr ett 1997). The 
flux of imag es A and B in the ra dio at 5 GHz was de- 
termined by IGarrett et al.l l|1994fl to be 15.8 mJy and 
17.2 mJy respectively. 

This is the highest redshift system in our sample, and 
is also interesting since one of the data points lands 
in the range where the 2175 A bump in the Galac- 
tic extinction law should be (see Figure I21|l . How- 
ever, the extinction signal is very weak with AiV) = 
0.1 ± 0.1, 0.11 ± 0.09, 0.28 ± 0.03 at xl = 1-9, 1.5, 1.6 
for the Galactic extinction law, power law and linear 
law respectively (see Table ^] for the parameters of the 
fits). When we fix Am we find somewhat higher extinc- 
tion of A{V) = 0.20 ± 0.03, 0.18 ± 0.03, 0.23 ± 0.01 at 
= 1.7, 1.4, 1.6. In both cases a power law or a linear 
law is marginally preferred to a Galactic extinction law. 
We also analyze the data with respect to a possible mi- 
crolensing signal but only find a weak microlensing signal 
(see Table [TBI and Figure [2J- Finally we note that our 
results for the Galactic extinction law are consistent with 
the results offFalco et al. (1999). 

4.2. The full sample 

In this section we study the properties of the sample 
as a whole. We look for correlations between various pa- 
rameters and, in particular, search for any dependence 
on the redshift or the morphology of the galaxies. Fur- 
thermore, we discuss the low Ry values found in SN la 
studies and the possible complementarity of lensing ex- 
tinction curve studies. We study, on the one hand, a 
'golden sample' and, on the other hand, we analyze the 
full sample. The 'golden sample' is defined to include 
the image pair with the strongest differential extinction 
for each lens. In addition Q0142— 100 is excluded from 
the 'golden sample' (see H4.1.5|l . The 'golden sample' 
therefore consists of eight pairs of images, of which seven 
have strong enough extinction to analyse the extinction 
curve. If not otherwise stated, the results apply to the 
full sample. 

4.2.1. A{V) as a function of distance from center of the 
lensing galaxy 

To study the distribution of A{V) as a function of dis- 
tance from the lens galaxy, we analyze the sample using 
two methods. First we plot, in Figure !^ the differential 
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TABLE 13 

Extinction curve fit results for MG0414+0534: Al 



Ima 


ges 


Extinction 


A{V) 






a 


Am 


s 


xi 


Al- 


B 


Galactic 


0.07 ± 0.02 





4±0.1 




-0.99 ± 0.03 




1.5 ±0.5 


Al- 


B 


Galactic 


0.15 ± 0.03 





9 ±0.2 




-1.07 




1.6 ±0.5 


Al- 


C 


Galactic 


0.29 ± 0.04 


1 


5 ± 0.2 




-2.00 ± 0.04 




1.7±0.5 


Al- 


C 


Galactic 


0.27 ±0.04 


1 


4 ±0.2 




-2.0 




1.5 ±0.4 


Al- 


B 


Power law 


0.10 ±0.05 






3.1 ±0.6 


-1.04 ±0.05 




0.9 ±0.4 


Al- 


B 


Power law 


0.13 ± 0.03 






2.8 ±0.4 


-1.07 




1.0 ±0.4 


Al- 


C 


Power law 


0.15 ± 0.07 






2.7 ±0.6 


-1.93 ± 0.05 




1.0 ±0.4 


Al- 


C 


Power law 


0.24 ± 0.04 






2.1 ±0.3 


-2.0 




1.1 ±0.4 


Al- 


B 


Linear law 


0.53 ± 0.04 






1.0 


-1.37 ±0.05 




2.2 ±0.5 


Al- 


B 


Linear law 


0.31 ± 0.02 






1.0 


-1.07 




3.4 ±0.4 


Al- 


B 


Linear law 


0.49 ± 0.03 






1.0 


-1.07 


0.28 ±0.05 


2.3 ±0.4 


Al- 


C 


Linear law 


0.61 ± 0.04 






1.0 


-2.27 ±0.05 




1.9 ±0.5 


Al- 


C 


Linear law 


0.42 ± 0.02 






1.0 


-2.0 




2.9 ±0.4 


Al- 


c 


Linear law 


0.58 ± 0.03 






1.0 


-2.0 


0.26 ±0.05 


1.9 ±0.5 



Note. — The extinction of imago Al compared to reference images B and C. Numbers quoted in 
italics were fixed in the fitting procedure. 



TABLE 14 

Extinction curve fit results for MG0414+0534: A2 



Images 


Extinction 




A{V) 


Rv 


a 


Am 


s 




A2-B 


Galactic 





87 ±0.05 


2.7 ±0.2 




-1.08 ±0.04 




1.8 ±0.5 


A2-B 


Galactic 





69 ± 0.03 


2.2 ±0.2 




-0.93 




1.8 ±0.5 


A2-C 


Galactic 





91 ± 0.04 


2.6 ±0.1 




-1.96 ±0.04 




1.8 ±0.5 


A2-C 


Galactic 





81 ± 0.04 


2.3 ±0.2 




-1.89 




1.6 ±0.5 


A2-B 


Power law 




1.6 ± 1.2 




0.7 ±0.2 


-1.8± 1.2 




1.7 ±0.5 


A2-B 


Power law 





65 ± 0.03 




1.49 ±0.07 


-0.9.3 




2.3 ±0.5 


A2-C 


Power law 




1.3 ±0.5 




0.9 ±0.2 


-2.4 ±0.5 




1.8 ±0.5 


A2-C 


Power law 





75 ± 0.03 




1.42 ±0.07 


-1.89 




1.9 ±0.5 


A2-B 


Linear law 


1 


11 ±0.04 




1.0 


-1.34 ±0.05 




1.7 ±0.5 


A2-B 


Linear law 





81 ±0.02 




1.0 


-0.9.3 




4.0 ±0.4 


A2-B 


Linear law 


1 


07 ± 0.04 




1.0 


-0.93 


0.39 ±0.05 


1.6 ±0.5 


A2-C 


Linear law 


1 


16 ± 0.04 




1.0 


-2.25 ±0.05 




1.6 ±0.4 


A2-C 


Linear law 





91 ±0.02 




1.0 


-1.89 




3.5 ±0.4 


A2-C 


Linear law 


1 


13 ±0.04 




1.0 


-1.89 


0.35 ±0.05 


1.5 ±0.4 


Note. 


— The extinction of imaj 


;e A2 compared to reference images B and C. Numbers quoted in 



italics were fixed in the fitting procedure. 



TABLE 15 
MG0414+0534: The extinction 

PROPERTIES OF A2 — Al 





Am 


A(V) 


Rv 


Fit 


Free 


0.61 ± 0.11 


3.8 ±0.7 


Fit 


Fixed 


0.53 ± 0.03 


3.5 ±0.4 


C 


Free 


0.62 ± 0.06 


4.0 ±0.9 


C 


Fixed 


0.54 ±0.06 


3.4 ± 1.0 


B 


Free 


0.80 ±0.05 


5.4 ± 2.5 


B 


Free 


0.54 ± 0.04 


3.7 ± 1.3 


Note. — 


The first two 


lines give 



the results from the Galactic extinc- 
tion law fit to the data. The last four 
lines give the exti nct ion properties cal- 
culated from eq. 1141 using the proper- 
ties of A2 and Al co mpa red to images 
B and C from Tables [T51 and [HI 



A{V) of the image pairs (from Table [T7|l as a function of 
the ratio of the distances from the center of the galaxy 
(from Table El). We assign a negative value to the A{V), 
in those cases where the more distant image is the more 
strongly extinguished one. One can see, that when the 



ratio is small, the image which is nearer the center of 
the galaxy is the more extinguished one. However, when 
the ratio approaches one, the A{V) becomes more evenly 
scattered around zero. 

For the second method, we assume that the image with 
the weakest extinction signal is indeed non-extinguished. 
We define an absolute A{V) for the other images, by 
taking the differential extinction compared to this refer- 
ence image, which we plot as a function of distance from 
the center of the galaxy, scaled by the lens galaxy scale 
length^ (see Figure ESI and Table [21). From the plot we 
can see that most A{V) values lie in the range of 0-0.5 
for distances smaller than around four scale lengths, but 
drop for more distant images. 

Both of these results are consistent with the expec- 
tation that the more distant image is on average more 
likely to pass outside the galaxy and thus not be affected 
by extinction. When the distances become similar, sec- 
ondary effects due to the non-symmetric shape of the 
lens start becoming important, creating a scatter in the 
A{V) vs. distance plots. This is in particular the case for 

^ The scale length is ta ken to be the effec tive radius of a de 
Vaucouleurs profile fit from lR.usin et all 120031^ . 
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TABLE 16 

Extinction curve fit results for MG2016+112 



Extinction 


A(V) 


Ry 


a 


Am 


s 




Galactic 


0.20 ±0.03 


3.0 ±0.5 




-0.092 




1.7 ±0.4 


Galactic 


0.1 ±0.1 


1.8 ± 1.0 




-0.01 ± 0.10 




1.9 ± 0.5 


Galactic 


0.11 ±0.09 


1.8 ± 1.3 




-0.092 


-0.1 ± 0.1 


1.8 ±0.5 


Power law 


0.18 ±0.03 




1.4 ±0.2 


-0.092 




1.4 ±0.4 


Power law 


0.11 ±0.09 




1.8 ±0.6 


-0.02 ± 0.10 




1.5 ± 0.5 


Power law 


0.12 ±0.09 




1.7 ±0.5 


-0.092 


-0.06 ±0.11 


1.5 ± 0.5 


Linear law 


0.23 ±0.01 




1.0 


-0.092 




1.6 ± 0.4 


Linear law 


0.28 ±0.03 




1.0 


-0.18 ± 0.06 




1.6 ± 0.5 


Linear law 


0.27 ±0.03 




1.0 


-0.092 


0.08 ±0.06 


1.6 ±0.5 



Note. — The extinction of image B compared to reference image A. Numbers quoted in 
italics were fixed in the fitting procedure. 
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Fig. 21. — MG2016+112: The upper panel shows the data points 
and the best fit extinction curves as given in Tableau The lower 
panel shows the original data points and their shift due to a mi- 
crolensing signal. The paramete rs o f the fits can be seen in Ta- 
ble ll6l See the caption of Figure ITTl for annotation overview. 




0.0 0.2 



Fig. 22. — The differential A(V') for a pair of images vs. the 
ratio of the distances from the center of the lensing galaxy. The 
differential A{V) is defined as negative if the image closer to the 
galaxy is less extinguished. The figure shows that images closer 
to the galaxy tend to be the more extinguished but that when the 
ratio approaches 1 the scatter increases. 
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Fig. 23. — AiV) as a function of the distance of the image relative 
to the scale radius of the lensing galaxy. We assign an absolute 
A(V) to the images by assuming that the least extinguished image 
for each system is non-extinguished. The quads are symbolized 
as circles and the four doubles are symbolized as triangles (up 
and down facing), a diamond and a box. The non-extinguished 
reference images are marked on the plot by open symbols. The 
error bars for dg/ds are smaller than the plotted symbols. We see 
that the A{V) values mostly lie in the range of — 0.5 for dg/ds ^ 4 
and drop for higher values. 



the quads where the distances tend to be similar. 

4.2.2. The different extinction laws 

We investigate whether our sample shows a preference 
for one type of extinction law to another and whether 
the type of extinction depends on the galaxy type. We 
also study the correlation between the parameters of the 
different fits. 

We find that when Am is allowed to vary, our sample 
does not show a preference for one extinction law over 
the other (the mean of the is — l-^, 1-6, 1.8 for the 
Galactic extinction law, power law and linear law) al- 
though individual systems can show a strong preference. 
If we alternatively look at the fits where Am was fixed we 
see that the power law and Galactic extinction law are 
preferred over the linear law in the sample as a whole 
(with xl = 2.7,2.1,4.3) but again individual systems can 
show different behaviors. 

There are three late type galaxies in our sample. One 
(HE0512— 3329) shows a clear preference for an SMC lin- 
ear law extinction, one (Bl 152-1-199) shows a preference 
for a Galactic extinction law and the third (Q2237+G3G) 
gives equally good fits to all the extinction laws (which 
is expected due to its low redshift, see S I2.6.1|l . For the 
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1.5 
A(V) 



2.0 2.5 3.0 



Fig. 24. — The power index a vs. Ry ■ The data points consist of 
all image pairs where the extinction curves were analyzed (from fits 
with Am free). The plot shows a clear correlation between a and 
Ry , with lower a giving higher Ry , which is consistent with lower 
Ry giving a steeper rise into the UV. The point corresponding to 
a = 1.0 (SMC type extinct ion) and Ry = 2.9 (the mean Ry for 
the SMC, as determined bv lPeiin-9921) . is marked by a diamond. 



early type galaxies there is also no clear preference for 
one type of extinction law. Three systems (PG1115+080, 
Q0142-100, RXJ0911+0551) show no preference for one 
extinction law over the other, one (MG0414+0534) fa- 
vors a power law with power index a — 2 — 3 for one of 
the images (which may be affected by extinction along 
both lines of sight) and one (MG0216-I-112) shows a weak 
preference for a power or a linear law over the Galac- 
tic extinction law. We therefore conclude that there is 
no evidence for a correlation between galaxy type and 
type of extinction in our sample. 

When confining the analysis to a Galactic extinction, 
we find that the mean Ry (for the 'golden sample') of the 
late type galaxies (i?^**= = 2.3 ± 0.5) is marginally lower 

than that of the early type galaxies {R'^^''^ = 3.2 ± 0.6), 
however they are consistent within the error bars and the 
difference may be due to low number statistics. This is 
further discussed in H4.2.4I 

We also study the correlation of the parameters Ry 
and a for each system to demonstrate the consistency 
of the two approaches. As expected, we find that there 
is a strong correlation with larger Ry giving smaller a, 
as seen in Figure]^ consistent with smaller Ry giving 
steeper rise into the UV. The exact relationship between 
Ry and a can be derived by solving equations (O and |(7J) 
and is wavelength dependent. A first order linear fit to 
our data gives a = (2.5 ± 0.2) - (0.51 ± 0.09)i?y which is 
an applicable approximation within the wavelength range 
of our data. In addition we check whether the strength 
of the extinction is correlated to Ry, but we find no 
evidence for such a correlation (see Figure [23). 

Finally we study the correlation between the values of 
A(V) which were found corresponding to the different ex- 
tinction laws. The results can be seen in Figure EEl We 
see that when Am is free, the power law favors higher 
A{V) than the other two extinction laws (see left bottom 
and top panels in Fig. l^ . When Am is fixed, the corre- 
spondence becomes much better with the power law giv- 
ing marginally lower values. The agreement between the 
A{V) values derived for the Galactic extinction law and 
the linear law are in general good, regardless of whether 
wc keep Am fixed or free (with the exception that for 
A{V) > 1 the linear law gives higher results when Am is 



Fig. 25. — The figure shows Ry vs. A{V). The data points 
consist of all image pairs where the extinction curves were analyzed 
(filled circles from fits with Am free and filled boxes from fits with 
Am fixed). The figure shows, as expected, no correlation between 
Ry and the amount of extinction. 




Fig. 26. — The graphs show the correlation between the val- 
ues of A{y) derived for the different fits. The left column shows 
the distribution for the fits where Am is a fitted parameter and 
the right column gives the corresponding distribution when Am is 
fixed. The dashed line corresponds to a; = j/, which would corre- 
spond to perfect agreement in Aiy) between the fits, and is plotted 
for reference. 



free) . 

Figure EHl also clearly demonstrates that the A{y ) val- 
ues become much better constrained when Am is fixed 
in the fitting. This suggests that it would be valuable for 
a future extinction survey, to do a simultaneous radio 
survey for the systems, in order to constrain the intrinsic 
magnitude ratio of the images. 

4.2.3. Evolution with redshift 

Next, we investigate the behavior of our sample as a 
function of redshift. The plots of Ry and a as a func- 
tion of redshift can be seen in the upper two panels of 
Figure [771 Wc do not see any strong correlation, in ei- 
ther the full nor the 'golden sample', although the lower 
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Fig. 27. — The top panel shows Rv as a function of z. The 
middle panels shows a vs. z. The bottom panels shows A(V) (as 
given by a Galactic extinction law fit) vs. z. On all panels, trian- 
gles denote late type galaxies and boxes denote early type galaxies 
where the values are taken from the fits with Am kept free. Circles 
denote the corresponding fits where Am was fixe d. Filled symbols 
correspond to a 'golden' sample as defined in il4.2l We see no strong 
evolution with z but lower values of Ry seem to appear at higher 
z. 



values of Rv and the higher values of a seem to appear 
at higher z. Our results do not confirm evolution of Ry 
with redshift, with lower Ry at higher z, as suggested by 
lOstman et all l)200()() . but do not exclude such evolution 
either. We stress that a larger sample would be needed 
to make any conclusive claims. 

We also investigate whether there might be an evo- 
lution in the amount of dust extinction with redshift. 
Again, we use two samples, the 'golden sample' as defined 
above, and a sample consisting of the highest differential 
extinction deduced for each image. The resulting plot 
can be seen in the bottom panel of Figure [23 and does 
not show any correlation between A{V) and z. 

4.2.4. Low Rv values and Type la SNe 

Recent studies of SNe la have suggested that Ry 
values for SN host s (which are mostly late type, see 
[Sullivan et al."2003') could be lower than those f or the 
Milky Way (see e.g., Riess et al. 1996b; Krisciunas et al.l 
UpOO; Wane 2006) suggesting that SN hosts a re systemat- 
ically different from the Milky Way. iWand l|2005^) how- 
ever suggest that the reason for the low values of Ry 
might be due to circumstellar dust around the SNe them- 
selves. The presence of such dust would cause inaccurate 
estimates of the dust extinction of the host galaxies of 
the SNe la. 

Lensing studies have also seen more extr eme Rv val- 
ues than thos e in the Milk y Way (see e.g.. iMotta et alJ 
H)02; Wucknit zerani200.lD but this has been criticized 
as possibly being due to extinction along both lines of 
sight. However, by choosing systems where the extinc- 
tion of the measured image dominates the extinction of 



the reference image, this effect can be avoided (see dis- 
cjussion in § I2.3|l . In our sample the mean Ry value is 
Ry = 2.8±0.3 (with RMS scatter of 1.2) for the full sam- 
ple, and Ry = 2.8 ± 0.3 (with RMS scatter of 1.1) for 
the 'golden sample'. These values are marginally lower 
than, but consistent with, the Milky Way mean value 
of Ry =3.1. If we look at Ry for the late and early 
type galaxies separately, we find R'^^'^^ = 3.2 ± 0.6 and 
j^iate _ 2 3 -|- Q 5 fQj. 'golden' sample. These values 
are consistent with each other within the quoted error 
bars, but it is interesting that the late type galaxies have 
a lower mean Ry, in agreement with SN la studies. A 
larger sample would be needed to determine whether this 
is a real trend, or due to low number statistics. 

The mean extinction in our 'golden sample' is A{V) — 
0.56 ± 0.04 (with RMS scatter of 0.80). This gives a 
lower limit to the mean absolute extinction, as the mean 
A{V) value is lowered if the reference image is also ex- 
tinguished. If we remove the highly extinguished system 
B1152-I-199 from our sample, the mean of the 'golden 
sample' reduces to A{V) = 0.29 ± 0.05 (with RMS scat- 
ter of 0.29). If we instead take the highest differential ex- 
tinction for each image for the full sample into account 
we get A{V) = 0.33 ± 0.03 (with RMS scatter of 0.56 
(or A{V) = 0.21 ± 0.03 without Bl 152+199). AU these 
values are high enough to cause systematic effects in the 
calibration of SNe la. 

A lower mean Ry from a lensing study would 
strengthen the results of low Ry values from SN studies 
applying to the interstellar medium. A lower real Ry 
value than the assumed one would lead to an overesti- 
mation of A{V) given a measurement of E{B — V) (as 
A{V) = RvE{B-V)). It is therefore important that the 
extinction properties of higher redshift galaxies, and SN 
hosts in particular, be further investigated as assuming a 
mean Galactic extinction with Ry = 3.1 in the analysis 
of SNe la could affect the cosmological results. 

Lensing galaxies and SN la hosts are distributed over 
a similar redshift range (from z = to 2: 1) and con- 
sist of both early type and late type galaxies. The ma- 
jori ty of lensing galaxies are massive early type galax- 
ies HKochanek et a1.ll2000l) whereas the SN la hosts are 
mostlv late tvpe l|Sullivan et al.l 12003(1 . It would how- 
ever be possible to select sub-samples of either group 
which would have the same morphology distribution as 
the other. Therefore, studies of the extinction properties 
of lensing galaxies can complement future dark energy 
SN la surveys, providing an independent measurement 
of the extinction properties of the SN la type hosts. 

5. SUMMARY 

We have presented an imaging survey of the extinction 
properties of 10 lensing galaxies using multiply imaged 
quasars observed with the ESO VLT in the optical and 
the NIR. We have made a dedicated effort to reduce the 
number of unknowns and effects which can mimic extinc- 
tion. We have explored, analytically and in simulations, 
the effects of extinction along both sight lines. We find 
that it is not crucial for the reference image to have zero 
extinction, as long as its extinction is small compared to 
the other image. We also study the effects of achromatic 
microlensing and find that to account for such an effect 
in photometric data, it is crucial to have constraints on 
the intrinsic magnitude difference of the images. 
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We were able to study the extinction of 9 out of 10 
of the systems in the siirvey, the last one had to be dis- 
carded due to contamination by the lensing galaxy. Out 
of the 9 systems, 8 have a two sigma extinction signal for 
at least one image pair, which was our limit for doing fur- 
ther extinction curve analysis. However, we suspect that 
one of those is also contaminated by the lensing galaxy 
(Q0142— 100) and exclude it from our 'golden sample'. 
The mean extinction for the 'golden sample' is A{V) = 
0.56 ± 0.09, using Galactic extinction law parametriza- 
tion^ and the mean Ry is Ry = 2.8 ± 0.4 (compared 
to Rv = 2.8 ± 0.3 for the full sample), which is con- 
sistent with the mean i?y — 3.1 found for the Galaxy. 
The systems show various extinction properties. There 
is no strong evidence for a correlation between morphol- 
ogy and extinction properties. As our sample covers a 
broad range in redshifts (z = 0.04 1.01) we have also 
looked for evolution with redshift. However, our results 
neither confirm nor refute evolution of extinction param- 
eters with redshift and we stress that a larger sample 
would be needed to make any conclusive claims. 

Finally we wish to point out that large studies of grav- 
itationally lensed quasars are ideal to study the possi- 



ble evolution of extinction properties as they are spread 
over a redshift range from z — to z 1. Furthermore, 
the quasars do not affect the environment of the galaxy 
we wish to study as is the case in SN la studies. For 
further improvements, however, higher resolution and 
deeper imaging for a larger sample would be required 
making a dedicated study with space based telescopes 
of considerable interest. A simultaneous radio survey, in 
order to constrain the intrinsic ratio of the images, would 
also further improve the results. Such a study could com- 
plement future dark energy SN la surveys, providing an 
independent measurement for the extinction properties 
of SN la type host galaxies. 
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Extinction curves of lensing galaxies 



TABLE 17 

Overview of extinction properties for the 

Galactic extinction law fit. 



Lens 


Lens redshift 


Image pair 


AiV) 


Ry 




0.04 


B— A 


\J,\J'-t _1_ \J.\J'-t 








C— A 


n 9Q _|_ n riR 


2.9 ± 1.4 






D— A 




3.1 ± 1.5 






C— B 


0.2 it 0.1 








D— B 


n 97 -1- n 1 

U.Z i ^ U- -LO 








D— C 


n no _|_ n f)"^ 
u.uo ni u.uo 






0.31 


A2— Al 


1 9 + n 0(S 

U. J-ii _1_ U.UU 


3.3 ± 1.9 






B— Al 


n + n 04 

U.UU _1_ U.Urt 








C— Al 


n -1- n (T^ 

u.uo ^ U.UU 








A2— B 


n fifi -1- n 0*^ 
U.UU ^ u.uo 








A2— C 


n 07 -1- n'^ 

U.U / ^ U.UU 








B— C 


n HQ _|_ n nc 
u.uo _i_ u.uo 




R14224-2'i1 


0.34 


B— A 


OS + 20 

u.uo _1_ U.tiU 








A— C 


0.1 di 0.3 








D— A 










B— C 


0.1 di 0.3 








D— B 










D— C 






R1 1 '^9-1-1 QQ 


0.44 


B— A 


9 4Q _|_ n fiQ 
^.lo ^ u.uy 


2.1 ±0.1 


001 49— 1 on 


0.49 


B— A 


n 40 _|_ n HQ 
u.iu ^ u.uo 


4.7 ± 0.7 




0.60 


B— A 








0.77 


B— A 


U . OO _1_ U.UU 


4.9 ± 0.6 
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n on _|_ n 1 fi 

U . .^U ^ U. -LU 


3.3 ± 1.5 






D-A 


0.12 ±0.06 








B-C 


0.02 ± 0.06 








B-D 


0.20 ±0.08 


4.1 ± 1.1 






C-D 


0.09 ± 0.08 


2.2 ± 1.5 


HE0512-3329 


0.93 


A-B 


0.14 ±0.04 


1.7 ±0.4 


MG0414+0534 


0.96 


A2-A1 


0.6 ±0.1 


3.8 ±0.7 






Al-B 


0.07 ±0.02 


0.4 ±0.1 






Al-C 


0.29 ± 0.04 


1.5 ±0.2 






A2-B 


0.87 ± 0.05 


2.7 ±0.2 






A2-C 


0.91 ±0.04 


2.6 ±0.1 
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0.01 ±0.06 




MG2016-1-112 


1.01 


B-A 


0.1 ± 0.1 


1.8 ± 1.0 



Note. — Table of extinction properties for the 10 lensing systems. 
The systems are ordered according to increasing redshift. The fit results 
are from the Galactic extinction law with Am free. 



